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■  It  is  common  practice  to  use  computer-generated  noise  contours  or  noise  zone  maps  to 
assess  noise  impact  and  perform  noise  related  land-use  planning.  However,  developers 
and  other  interested  parties  often  question  the  accuracy  of  computer  simulations  and 
suggest  direct  measurement  to  “verify”  the  computer  predictions.  This  report  quantifies 
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STRATEGIES  FOR  AND  VALIDITY  OF 
NOISE  MONITORING  IN  THE  VICINITY 
OF  CIVILIAN  AIRFIELDS  AND 
ARMY  INSTALLATIONS 


1  INTRODUCTION 


Background 

It  is  common  practice  to  use  computer-generated 
noise  contours  or  noise  zone  maps  to  assess  noise 
impact  and  perform  noise-related  land-use  planning. 
In  the  United  States,  noise  zone  maps  are  usually 
expressed  in  terms  of  the  day/night  average  sound 
level  (DNL)  descriptor.1  Most  noise  zone  maps  are 
created  by  computer  simulation  programs  like  the 
Federal  Aviation  Administration’s  (FAA’s)  Integrated 
Noise  Model  (INM),  the  Air  Force’s  NOISEMAP, 
or  the  Army’s  BNOISE.3 

When  noise  zone  maps  are  used  for  noise  assess¬ 
ment  and  especially  land-use  planning,  developers 
and  other  interested  parties  often  question  the  accu¬ 
racy  of  the  computer  simulations  and  suggest  direct 
measurement  to  “verify”  the  computer  predictions. 
It  is  naturally  assumed  that  direct  measurement  must 
be  more  accurate  than  computer  simulation. 

Purpose 

The  purpose  of  this  report  is  to  quantify  the  tempo¬ 
ral  sampling  requirements  for  and  the  accuracy  and  the 
ability  of  directly  measured  sampled  data  to  estimate 
the  true  yearly  DNL. 

Approach 

The  Army’s  main  concern  is  the  blast  noise  created 
by  such  operations  as  armor,  artillery,  and  demolition. 
However  most  existing  monitor  data  were  gathered 
near  major  metropolitan  airports.  Thus,  the  analysis  of 


1  Guidelines  for  Considering  Noise  in  Land  Use  Planning 
and  Control  (federal  Interagency  Committee  on  Urban  Noise, 
June  1980). 

* Integrated  Noise  Model  (INM)  (Department  of  Trans¬ 
portation.  federal  Aviation  Administration);  Community 
Noise  Exposure  Resulting  From  Aircraft  Operations:  Com¬ 
puter  Program  Description.  AMRL-TR-73-106  (Department 
of  the  Air  force,  November  1974);  and  V.  Pawlowika  and  L. 
Little,  The  Blase  Noise  Prediction  Program:  User  Reference 
Manual,  Technical  Report  N-75/ADA074050  (U.S.  Army 
Construction  Engineering  Research  Laboratory  (OERL] , 
August  1979) 


monitoring  accuracy  with  specific  reference  to  blast 
noise  proceeds  along  two  paiallcl  paths.  For  the  first 
path,  the  results  of  limited  bias!  noise  monitoring  neai 
two  Army  installations  arc  compared  with  computci 
simulation  results  and  the  results  of  altitudinal  surveys 
of  the  community  response  to  blast  noise  in  the  areas 
near  these  same  installations.  For  the  second  path, 
metropolitan  airport  data  are  used.  These  data  exist 
for  a  number  of  airports  where  continuous  daily 
monitoring  was  performed  for  1  year  or  more  at  several 
sites.3  Studies  and  analyses  are  performed  on  these 
data  to  show,  quantitatively,  the  accuracy  that  dif¬ 
ferent  sampling  strategies  would  have  achieved. 

Attitudinal  surveys  are  used  to  gauge  or  quantify 
the  community  response  to  some  stimulus,  such  as 
noise.  During  the  past  30  years,  many  attitudinal 
surveys  have  been  conducted  worldwide  to  better 
understand  and  assess  human  and  community  response 
to  noise.  These  studies,  which  concentrate  mainly  on 
automobile  and  truck  traffic  and  rail  and  fixed-wing 
aircraft  noise,  have  resulted  in  a  proliferation  of  noise 
assessment  models  or  descriptors.  In  general,  these 
descriptors,  in  one  fashion  or  another,  take  into 
account  the  following: 

1 .  Sound  level  of  the  noise  events. 

2.  The  frequency  of  the  occurrence  of  the  noise 
events. 

3.  The  time  of  day  at  which  the  noise  occurs. 

In  general,  one  major  purpose  of  the  attitudinal 
survey  is  to  develop  a  highly  correlated  functional 
relation  between  some  measure  of  community  annoy¬ 
ance  (the  dependent  variable)  with  one  or  another  of 
these  noise  descriptors.  Over  the  past  few  years,  the 
scientific  community  has  generally  settled  on  the  use 
of  “high  annoyance”  as  a  measure  of  the  community 
response,  and  the  use  of  the  day/night  average  sound 
level  (DNL)  as  the  noise  descriptor.  High  annoyance  is 
defined  to  be  those  respondents  in  an  attitudinal  survey 


’Richard  E.  DeVor,  el  at..  “Development  of  Temporal 
Sampling  Strategies  for  Monitoring  Noise,  ”  Journal  of  the 
Acoustical  Society  of  America,  Volume  66,  No.  3  (September 
1979),  pp  763-771;  Paul  D.  Schomer  and  Richard  E.  DeVor, 
"Temporal  Sampling  Requirements  for  Estimation  of  Long- 
Term  Average  Sound  Levels  in  the  Vicinity  of  Airports," 
Journal  of  the  Acoustical  Society  of  America.  Volume  69, 
No.  3  (March  1981),  pp  713-719;  and  Paul  D.  Schomer.  et 
al.,  "Sampling  Strategies  for  Monitoring  Noise  in  the  Vicinity 
of  Airports."  Journal  of  the  Acoustical  Society  of  America. 
Volume  73.  No.  6  (June  1983),  pp  2041-2050.  (Note:  These 
papers  arc  included  as  Appendices  A.  B.  and  ('  of  this  report.) 


who  rate  themselves  in  the  top  two  annoyance  cate¬ 
gories  on  a  five-point  adjectival  scale  when  responding 
as  to  their  overall  annoyance  to  the  noise  environment. 
The  five  adjectival  ratings  are  (1)  extremely  annoyed, 
(2)  very  much  annoyed,  (3)  moderately  annoyed, 
(4)  a  little  annoyed,  (5)  and  not  at  all  annoyed.  The 
independent  variable  (the  noise  environment  described 
by  DNL)  is  normally  either  predicted  by  computer 
simulation  or  directly  measured.  Frequently,  direct 
measurements  are  used  to  spot  check  computer  simu¬ 
lations.  As  indicated  above,  one  purpose  of  this  report 
is  to  investigate  the  applicability  of  direct  measurement 
to  check  computer  simulation. 

For  computer  simulation,  every  noise  source  is 
tabulated  and  evaluated:  every  weapon  firing,  every 
target,  every  shell,  etc.  The  directivity  of  each  weapon 
for  each  firing  is  factored  and  a  statistical  distribution 
of  the  received  amplitudes  as  a  function  of  distance 
which  reflects  variable  weather  conditions  is  developed 
for  each  event.  The  received  noise  is  added  at  a  grid  of 
rectangular  points  covering  the  installation  and  the 
surrounding  area.  The  incremental  noise  from  every 
source  is  added  at  each  of  these  grid  points.  When  the 
summation  process  is  completed,  equal  noise  contours 
are  developed. 

Impulse  noise -the  noise  generated  by  armor, 
artillery,  or  demolition-is  assessed  using  the  C- 
frequency  weighting.  Noise  events  that  generate  sound 
levels  which  fall  below  a  certain  threshold  are  discarded. 
Both  direct  measurement  and  computer  simulation 
must  take  into  account  this  threshold  and  the  C- 
weighting.  Clearly,  the  computer  simulation  includes 
only  blast  noise  produced  by  the  Army  installation. 
For  direct  measurement,  care  must  be  exercised  to 
ensure  that  the  noise  monitoring  results  include  only 
Army  impulse  noise  and  not  other  noise  such  as  wind¬ 
generated  noise,  helicopter  flybys,  diesel  trains  or 
other  nearby  sources  of  high-level,  C-weighted  noise. 

Figure  1  shows  data  the  U.S.  Army  Construction 
Engineering  Research  Laboratory  (C ERL)  has  gathered 
which  relate  high  annoyance  and  the  computer- 
predicted  C-weighted  DNL  in  the  vicinity  of  Fort 
Bragg  and  Fort  Lewis.  This  figure  also  indicates  the 
National  Academy  of  Science's  recommended  function 
for  relating  community  response  to  C-weighted  DNL. 


These  data  indicate  that  if  the  noise  environment  at  an 
installation  is  specified  by  computer-predicted  DNL. 
then  the  resulting  annoyance  predicted  by  the  National 
Academy  of  Science  curve  will  likely  underestimate  the 
true  annoyance  by  a  small  amount. 

In  analyzing  the  viability  of  using  directly  monitored 
impulse  noise  data,  this  report  looks  at  the  results  of 
the  altitudinal  surveys  and  the  community  response  in 
terms  of  high  annoyance  as  rated  against  measured 
DNL  levels  and  the  National  Academy  of  Science 
recommendations.  The  bottom  line,  in  terms  of  the 
Army’s  interest,  is  accurate  prediction  of  the  com¬ 
munity  response  to  the  noise  environment  around 
Army  installations.  When  measured  DNL  improves  the 
Army’s  ability  to  properly  predict  community  annoy¬ 
ance,  direct  measurements  should  be  used  to  augment 
computer  prediction.  However,  when  direct  measure¬ 
ment  correlates  less  well  than  computer  prediction 
with  community  response  and  depreciates  the  Army's 
prediction  capability,  greater  reliance  must  be  placed  on 
computer  prediction  and  less  on  direct  measurement. 

As  an  approach  to  quantifying  temporal  sampling 
requirements,  approximately  1  year  of  daily  DNL 
(or  CNEL)  data  were  obtained  from  several  major 
metropolitan  commercial  airports.  These  data  were 
gathered  by  the  airports  using  several  fixed  monitors 
at  various  locations  around  the  airports.  The  time- 
series  data  are  modeled  and  various  statistical  analyses 
are  performed  on  these  data  to  show,  quantitatively, 
the  accuracy  that  different  sampling  strategies  would 
have  achieved. 

Mode  of  Technology  Transfer 

This  report  develops  guidelines  as  to  when  direct 
measurement  can  be  used  to  predict  community 
response  in  the  vicinity  of  Army  installations  and  when 
only  computer  simulation  should  be  used.  Temporal 
sampling  strategies  are  developed  for  use  with  direct 
measurement.  The  results  of  this  report  will  be  used 
by  the  Army  Environmental  Office  and  the  Army 
Environmental  Hygiene  Agency  to  formulate  their 
strategy  for  direct  measurement  as  a  part  of  Instal¬ 
lation  Compatible  Use  Zone  (1CUZ)  program  mandated 
by  Army  Regulation  (AR)  200-1.  Environmental 
Protection  and  Enhancement  (Department  of  the 
Army,  15  June  1982). 
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Figure  1.  Fort  Bragg  and  Fort  Lewis  data  (based  on  computer  prediction  only). 


2  ARMY  INSTALLATION 
NOISE  MONITORING 


In  1978  and  1980,  CERL  administered  attitudinal 
surveys  of  the  community  response  to  blast  noise  in 
the  vicinity  of  Fort  Bragg,  NC,  and  Fort  Lewis,  WA.4 
At  the  same  time  the  surveys  were  administered, 
CERL  extensively  monitored  the  actual  blast  noise 
produced  by  installation  operations  and  prepared 
computer  predictions  for  blast  noise  in  the  areas  sur¬ 
rounding  the  installations.  At  both  installations,  the 
computer  predictions  correlated  well  with  the  measured 
community  response.  The  results  of  direct  noise 
monitoring  at  Fort  Bragg  generally  correlated  well 
with  the  computer  predictions,  except  for  areas  distant 
from  the  installation  and  to  the  west  or  southwest. 

The  Fort  Bragg  study  was  performed  in  1978,  and 
the  Fort  Lewis  study  was  performed  in  1980.  Thus, 
this  chapter  first  describes  the  Fort  Bragg  results,  then 
the  Fort  Lewis  results.  The  details  of  these  studies  are 
discussed  below;  but  based  on  the  lessons  learned  at 
Fort  Bragg,  the  attitudinal  survey  at  Fort  Lewis  was 
conducted  at  eight  clusters  around  eight  fixed  monitor¬ 
ing  sites.  This  design  was  chosen  to  minimize  move¬ 
ment  of  the  monitors.  To  increase  reliability,  all 
monitors  at  Fort  Lewis  were  powered  from  110-V 
sources  and  the  monitoring  was  performed  for  a 
6-month  period. 

Fort  Bragg  Monitoring  Results 

In  the  vicinity  of  Fort  Bragg,  24-hour  monitoring 
was  done  at  17  sites.  The  number  of  complete  24-hour 
days  of  monitoring  at  each  site  ranged  from  4  to  67, 
with  25  being  a  typical  value.  All  C-weighted  data 
recorded  by  CERL’s  monitoring  equipment  were 
extensively  tested  and  checked  to  eliminate  all  but 
blast  noise.  All  data  were  recorded  in  6-minute  blocks. 
To  reduce  the  effects  of  noise  generated  by  wind  at 
the  microphone,  the  monitors  were  turned  off  when 
the  wind  meter  indicated  the  winds  were  blowing  at 
more  than  about  18  km/hour.  Whenever  the  monitors 
went  above  the  preset  peak  level  threshold  of  105 
dB  (95  dB  at  night),  an  analog  tape  recorder  and  a 
special  digital  timer  were  turned  on.  If  the  wind 
threshold  signal  came  on  at  any  time  when  the  recorder 
and  timer  were  running,  the  data  in  that  6-minute 


‘Paul  D.  Schomer,  Community  Reaction  to  Impulse 
Noise :  Initial  Army  Survey,  Technical  Report  N-I00/ADA 
101674  (CERL,  June  1981);  and  Paul  D.  Schomer.  Com¬ 
munity  Response  to  Impulse  Noise.  A  10-  Year  Research 
Summary.  Technical  Report  N-167  (CERL.  November  1983) 


block  were  discarded.  If  the  threshold  was  exceeded 
for  more  than  2  seconds,  a  technician  would  listen  to 
the  analog  tape  to  see  if  the  signal  was  caused  by 
impulses  or  by  some  other  source,  such  as  an  aircraft 
or  a  helicopter.  If  the  technician  detected  any  other 
type  of  noise  source  on  the  analog  tape,  that  6-minute 
data  block  was  also  discarded.  Thus,  the  only  data 
included  were  those  for  which  (1)  the  wind  threshold 
was  not  triggered  and  (2)  no  other  source  could  be 
heard  or  the  event  was  less  than  2  seconds  long,  or 
both. 

Figure  2  shows  the  general  outline  of  be  Fort 
Bragg  study  area,  overlaid  with  comput*-  ed'eted 
C-weighted  DNL  (CDNL)  contours  for  the  <r  before 
CERL’s  study.  Also  shown  are  15  of  the  17  .  toring 

sites  (two  sites  near  the  airfields  measured  aircraft 
noise  and  are  not  shown).  The  figure  grou  .:  ’  areas 
by  their  geographic  area  and  noise  zone,  l  -  a  off- 
installation  areas  in  the  same  general  region  and  noise 
zone  are  grouped  separately.  Table  1  lists  the  computer- 
predicted  and  directly  measured  noise  levels  by  moni¬ 
toring  site. 

At  sites  1  and  2  the  measured  levels  were  much 
higher  than  predicted  because  units  assigned  firing 
points  within  1  km  of  these  monitors  actually  fired 
from  much  closer  than  1  km.  In  the  areas  to  the  east, 
the  monitored  results  ranged  from  1 1  dB  below'  to  3 
dB  above  prediction.  For  the  sites  at  which  the 
measured  levels  were  close  to  prediction  in  the  east, 
the  predominant  noise  all  came  in  one  to  several  days, 
each  day  characterized  by  a  period  of  high  noise 
caused  by  sound-focus  conditions.*  In  contrast, 
monitor  sites  5,6,  8,  and  9,  to  the  south  and  west  of 
the  study  area,  exhibited  no  such  focus  days.  As  a 
result,  Table  1  shows  a  much  larger  difference  between 
the  computer-predicted  and  the  measured  values  for 
those  locations. 

The  differences  between  prediction  and  measure¬ 
ment  seem  to  follow  a  trend.  Sites  1  and  2,  which  were 
very  close  to  firing  points,  had  measured  data  which 
were  well  above  prediction.  Sites  3  and  4,  which  were 
about  1  mile  (1.6  km)  from  the  nearest  firing  point, 
had  measured  data  which  were  2  to  4  dB  above  pre¬ 
diction.  Sites  to  the  east  (both  on  and  off  the  instal¬ 
lation)  had  measured  data  which  were  somewhat  below 

•The  velocity  of  sound  changes  with  altitude  primarily 
because  of  changes  in  wind  velocity  and  temperature  with 
altitude.  This  sound  velocity  profile  can  focus  sound  much  as 
a  lens  focuses  light.  The  result  is  the  possibility  of  very  loud 
sounds  focused  at  far  distances  (e  g.,  2  to  25  miles)  from  the 


source. 


or  at  the  predicted  level.  Sites  to  the  south  and  west 
had  measured  data  which  were  far  below  the  predicted 
levels.  The  very  high  levels  at  Sites  1  and  2  are  believed 
to  have  been  caused  by  Marine  units  which  fired  from 
other  than  the  locations  they  listed.  (Large  percentage 
errors  in  small  distances  and  firing  points  closer  than 
300  m  to  monitors  are  beyond  the  scope  of  CERL’s 
noise  contour  prediction  computer  program.) 

For  those  sites  where  the  measured  data  agree  with 
prediction,  most  of  the  sound  energy  comes  during 
1  to  2  hours  over  a  few  days  when  focus  conditions 
existed  that  would  cause  high  noise  levels.  During 
other  times,  the  monitors  measured  much  lower 
noise  levels.  These  results  are  in  accordance  with  the 
statistical  nature  of  sound  propagation  resulting  from 
the  extreme  variations  caused  by  weather  conditions. 
The  sites  to  the  west  and  southwest  measured  levels 
well  below  prediction.  At  the  time  of  the  Fort  Bragg 
study,  it  was  thought  that  this  occurred  because  the 
monitoring  was  not  performed  for  a  long  enough  time. 

Fort  Lewis  Monitoring  Results 

At  Fort  Lewis,  24-hour  monitoring  was  done  at 
eight  sites  for  6  months  before  the  attitudinal  survey 
was  administered  (11  January  to  30  June  1980).  In 
comparison,  the  monitoring  duration  at  Fort  Bragg 
was  about  5  months  at  17  sites,  most  of  which  alter¬ 
nated  weekly.  To  minimize  the  equipment  problems 
encountered  at  Fort  Bragg  associated  with  moving 
monitoiing  locations,  eight  fixed  locations  were  chosen 
for  the  Fort  Lewis  study.  All  of  the  Fort  Lewis  sites 
were  powered  by  I  10-V  lines  and  included  the  normal 
uninterruptible  power  supply  in  the  (  ERL  monitor. 
As  at  Fort  Bragg,  the  data  were  tested  to  ensure  that 
only  blast  noise  was  included  in  the  monitored  results. 
The  same  wind  speed  and  peak  amplitude  thresholds 
were  used  at  Fort  Lewis  as  at  Fort  Bragg.  This 
modified  equipment  setup  improved  the  overall  per¬ 
formance  of  the  study  by  increasing  the  number  of 
successfully  completed  monitoring  days  by  at  least 
100  percent  and  reducing  equipment  failures  by  at 
least  200  percent. 

Figure  3  shows  the  general  outline  of  the  Fort 
Lewis  study  area  overlaid  with  computer-predicted 
CDNL  contours  for  1  year  preceding  the  monitoring 
period.  The  figure  also  shows  the  eight  monitoring 
sites:  these  sites  were  chosen  to  include  a  range  of 
community  types  (i.e.,  small  town,  city,  .uburban, 
and  cantonment  area). 

Table  2  lists  the  computer-predicted  and  measured 
CDNLs  by  monitor  location.  Although  the  predicted 


and  monitored  levels  tend  to  correlate,  they  do  not 
agree.  With  one  exception,  the  measured  levels  fall 
far  below  prediction.  These  same  results  were  found 
to  the  south  and  west  at  Fort  Bragg.  Apparently,  the 
poor  Fort  Bragg  results  to  the  south  and  west  were  not 
due  merely  to  “not  measuring  long  enough,"  as  had 
been  assumed. 

Although  a  number  of  theories  can  be  advanced  to 
explain  these  poor  results,  none  by  itself  gives  a  satis¬ 
factory  answer.  These  theories  include: 

1.  Blast  noise  was  lost  by  deleting  data  when  high 
winds  occurred  (above  18  km/hour). 

2.  Sound  focusing  conditions  tailed  to  occur  at 
the  eight  Fort  Lewis  monitor  sites  during  the  measure¬ 
ment  period. 

3.  The  peak  threshold  (105  dB  for  daytime  and  95 
dB  for  night)  deleted  meaningful  data. 

4.  The  monitoring  equipment  is  fundamentally 
incorrect  in  its  operation. 

5.  The  monitored  results  are  correct  and  the  pre¬ 
diction  and  measures  used  to  assess  community  response 
are  incorrect. 


Table  2 

Measured  vs  Predicted  CDNL  at  Fort  Lewis 


Site* 

Computer- 
Predicted 
During  the 
Monitoring 
Period** 
(CDNL) 

Measured 

(CDNL) 

Difference 

Computer- 
Predicted 
for  the 
Entire  Year 
(CDNL) 

1 

44 

34 

+  10 

45 

2 

45 

30 

+  15 

47 

3 

59 

46 

+  13 

61 

4 

53 

46 

+7 

53 

5 

56 

43 

+  13 

56 

6 

53 

37 

+  16 

53 

7 

49 

53 

4 

50 

8 

54 

44 

+  10 

54 

*See  I  ifjurc  3  lor  site  locations. 

*  *  All  stations  ran  tor  the  entire  6  months  with  only  a  lew 
scattered  days  of  data  lost. 


s 

* 

t 


Figure  3.  Predicted  CDNL  contours,  monitor  sites  and  predominant  respondent  groups  in  the  Fort  Lewis  study  area 


Figure  4  arrays  the  at  tit  udinal  survey  results  in 
terms  of  community  response  (i.e.,  high  annoyance)  vs 
the  monitored  DNL.  This  figure  also  includes  the 
National  Academy  of  Science’s  recommended  function 
for  predicting  community  response  based  on  CDNL. 
Clearly,  if  a  measured  DNL  value  is  plugged  into  the 
recommended  National  Academy  of  Science  relation 
for  predicting  community  response,  then  community 
annoyance  is  greatly  underestimated  as  compared  with 
the  results  from  the  attitudinal  surveys.  In  the  discus¬ 
sion  which  follows,  this  result  (i.e.,  that  annoyance 
prediction  based  on  measured  levels  greatly  under¬ 
estimates  true  annoyance)  and  the  corresponding 
results  for  computer  simulation  and  the  data  in  Figure 
1  are  used  to  analyze  some  of  the  causes  to  explain  the 
discrepancy  between  results  obtained  by  computer 
simulation  and  by  direct  measurement. 


Fort  Bragg  and  Fort  Lewis  Results  -  Discussion 

1 .  Can  the  monitoring  be  correct  and  the  computer 
predictions  and  the  attitudinal  survey  results  be  in¬ 
correct?  The  following  facts  are  known:  training, 
including  artillery  fire,  mortars,  and  demolition, 
occurred  at  a  more  or  less  normal  rate  (perhaps  at  j 
somewhat  decreased  rate)  during  the  Fort  Lewis 
monitoring  period.  The  attitudinal  survey  results  in 
terms  of  community  annoyance  correlate  well  with  the 
predicted  noise  environment  and  are  in  general  agree¬ 
ment  with  the  previous  survey  results  and  prediction 
at  Fort  Bragg.  Thirty  percent  of  about  1500  respon¬ 
dents  interviewed  at  Fort  Lewis  report  hearing  blast 
noise  either  daily  or  several  times  per  week,  and  two- 
thirds  of  these  respondents  say  that  the  blast  noise 
is  much  louder  than  ordinary  conversation.  The  survey 
interviewers  also  occasionally  reported  that  they  heard 
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Figure  4.  Survey  data,  percent  highly  annoyed  (HA)  vs  measured  CDNL.  The  dashed  line  is  the  CHABA  recom¬ 
mended  relation. 


blasl  noise  when  I  hey  conducted  interviews;  sometimes  Tins  is  the  threshold  set  as  the  recommended  practice 

during  an  interview,  the  building  would  shake  from  of  the  National  Academy  of  Science  Committee  on 

blast  noise  excitation.  These  facts  seem  to  indicate  Hearing,  Biocoustics  and  Biomechanics  (CHABA). 

that  the  answer  cannot  lie  wholly  in  the  statement  Mathematically,  there  would  have  to  be  between  100 

that  no  noise  existed.  Rather,  it  appears  there  was  a  and  1000  of  these  "just  missed"  events  per  day  to 

failure  to  properly  record  some  of  the  blast  noise  data.  bring  the  measured  values  into  general  agreement  with 

the  predicted  levels.  So,  while  it  also  may  be  true  that 

2.  Does  the  noise  monitoring  equipment  properly  the  CHABA  procedure  should  not  incorporate  this 

measure  impulse  noise?  CERL  has  had  about  4  years  of  85-dB  threshold,  the  numbers  are  so  small  that  the 

experience  with  this  equipment.  It  has  been  tested  in  total  answer  cannot  lie  with  the  threshold.  In  the 

comparison  with  every  major  commercial  noise  case  of  Fort  Lewis,  there  are  not  enough  firings  to 

monitor  and  always  equals  or  betters  any  commercial  even  theoretically  explain  the  discrepancy  on  the 

equipment  in  terms  of  measuring  a  known  noise  basis  of  the  85-dB  threshold. 

source.  It  has  measured  blast  noise  at  Fort  Bragg  and 

Fort  Carson  using  the  same  general  techniques.  It  has  5.  Did  all  of  the  monitors  fail  to  he  located  at 

been  used  successfully  lor  Army  source  noise  measure-  focus  points ?  The  sound  levels  received  in  the  com¬ 
ments  and  for  many  Army  sound  propagation  measure-  mimily  vary  statistically.  Weather  conditions  can 

ments.  It  is  designed  and  operated  in  accordance  with  focus  high-amplitude  sound  at  distant  locations.  For 

all  applicable  American  National  Standards  Institute  example,  the  Army  environmental  Hygiene  Agency 

standards.  A  basic  flaw  in  the  equipment  or  general  recorded  an  instance  of  a  noise  complaint  at  one 

operation  does  not  appear  to  be  a  reasonable  explan-  installation  where  a  woman  reported  her  house  shaking 

ation  for  the  poor  results  obtained  from  these  studies.  from  artillery  noise.  Sometime  within  the  20  minutes 

it  took  someone  from  the  installation  to  go  out  to  the 

3.  Does  deletion  of  windy  data  delete  important  woman’s  house,  her  house  stopped  shaking.  But  the 

blast  noise  data?  The  general  monitoring  procedure  bam  a  few  hundred  feet  away  from  the  house,  had 

used  at  Fort  Bragg  and  Fort  Lewis  was  to  delete  all  begun  shaking.  The  point  of  this  anecdote  is  that  sound 

blocks  of  data  for  which  the  wind  exceeded  18  km/  focuses  can,  at  times,  be  very  localized  in  nature.  It 

hour.  This  was  done  because  wind  turbulence  gen-  may  be  tfl>,  none  of  (he  eight  For[  Uwis  monitors 

erated  at  the  microphone  can  cause  readings  which  regu,arly  received  these  sharply  focused  sounds 

will  appear  as  blast  noise.  The  known  facts  indicate  a]though  some  percentage  of  the  population  near  these 

the  following,  there  were  not  many  windy  days  during  monitors  did  receive  the  sharply  focused  sounds, 

the  monitoring  period  at  Fort  Lewis,  so  little  data  were  However,  this  explanation  does  not  seem  very  likely, 

deleted  because  of  the  wind  threshold.  On  the  other 

hand,  blast  propagation  measurements  conducted  by 

CBRL  at  Fort  Leonard  Wood  showed  that  the  greatest  Summary 

sound  propagation  occurred  on  1  day  during  the  4-  directly  measured  blast  noise  results  at  Fort 

week  test  period  when  a  wind  shear  occurred  at  about  Lewis,  and  to  some  extent  Fort  Bragg,  do  not  agree 

1000  ft  (305  m)  above  ground  level.5  Thus,  it  may  be  with  computer  prediction  or  with  the  community 

that  some  blast  noise  data  were  lost  because  of  the  resPonsc  in  terms  of  annoyance  or  reported  frequency 

wind  threshold.  However,  the  respondent-reported  of  hearin8  blast  noise-  Many  reasons  can  be  advanced 

frequency  of  hearing  blast  noise  substantially  exceeds  for  these  discrepancies,  but  none  appear  to  fully 

the  rate  at  which  windy  days  occurred  at  Fort  Lewis  e*Plain  them-  However,  three  of  the  potential  sources 

during  the  monitoring  period.  So  the  full  answer  of  discrepancy  can  be  mitigated: 

cannot  lie  in  deletion  of  “windy  data.” 

1.  The  wind  threshold  and  the  peak  amplitude 

4.  Can  the  104-dB  peak  threshold  be  deleting  blast  threshold  can  be  altered  if  windscreens  are  improved 

noise  data?  The  105-dB  peak  threshold  corresponds  to  and  if  a  multiple  microphone  technique  can  be  devel- 

a  blast  having  about  an  85-dB  sound  exposure  level.  oped  to  separate  wind-induced  noise  from  blast  noise. 

2.  With  windscreen  improvements  and  a  multiple 
microphone  array,  it  may  be  possible  to  set  the  peak 
threshold  at  95  or  100  rather  than  105  dB;  the  wind 
speed  deletion  level  could  then  be  set  at  perhaps 
30  instead  of  1 8  km/hour. 


’  P  D.  Schomer,  et  at..  The  Statistics  of  Amplitude  and 
Spectmm  of  Blasts  Propagated  in  the  Atmosphere.  Technical 
Report  N-13,  Volume  1  (ADA03347S)  and  Volume  2  (ADA 
03361)  «T:RL.  November  1976). 
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3.  The  only  answer  to  the  question,  "Are  we 
measuring  long  enough  to  obtain  a  good  statistical 
sample  of  what  the  weather  effects  are  likely  to  be  on 
the  sound  propagation?"  is  to  measure  for  a  very  long 
time,  and  clearly  this  can  be  done. 

The  general  results  of  the  Fort  Bragg  and  Fort 
Lewis  studies  seem  to  dictate  that  no  more  monitoring 
be  performed  for  noise  contour  verification  purposes 
(except  possibly  near  to  the  sources)  until  the  wind¬ 
screens  used  with  the  noise  monitoring  equipment  are 
improved  and  until  a  multiple  microphone  technique 
is  developed  which  can  better  separate  wind  effects 
from  true  blast  noise  data.  Also,  any  future  monitoring 
should  be  done  for  at  least  1  year  in  order  to  better 
account  for  the  extreme  variation  in  sound  propaga¬ 
tion  and  focus  location. 


SAMPLING  STRATEGIES  FOR 
3  MONITORING  NOISE  IN  THE 

VICINITY  OF  CIVILIAN  AIRPORTS 

Introduction 

When  one  can  monitor  blast  noise,  such  as  near  an 
installation,  then  a  temporal  sampling  strategy  must  be 
developed.  This  chapter  uses  commercial  airport  data 
to  develop  a  notion  about  temporal  sampling  require¬ 
ments  in  general.  The  next  chapter  relates  the  airport 
results  to  the  Army  situation. 

The  general  problem  underlying  temporal  sampling 
requirements  for  estimating  long-term  average  sound 
levels  at  civilian  airports  is  the  associated  statistical 
assessment  of  the  precision  of  the  estimates  of  mean 
sound  level.  With  only  a  few  exceptions,  most  tech¬ 
niques  in  use  today  for  sampling  community  noise 
call  for  sampling  over  relatively  short  periods  of  time, 
i.e„  from  a  few  minutes  to  perhaps  a  single  day. 
However,  the  time  varying  nature  of  noise  data  when 
viewed  as  a  time  series  (hourly  or  daily  averages) 
suggests  that  short-term  sampling  may  lead  to  serious 
inaccuracies  in  the  estimation  of  a  long-term  (yearly) 
mean  noise  level.  For  example,  the  24-hour  periodic 
pattern  in  hourly  mean  sound  level  may  vary  from 
about  40  to  85  dB.  The  Community  Noise  Equivalent 
Level  (CNEL)  and  DNL  noise  descriptors  both  com¬ 
monly  vary  from  45  to  80  dB.  These  wide  ranges  for 
sound  level,  together  with  the  fact  that  the  data  in 
general  exhibit  high  positive  autocorrelation  and  high 
coefficients  of  variation,  suggest  that  small  or  short 
sampling  periods  or  both  may  provide  imprecise  and 
inaccurate  mean  value  estimates. 


The  techniques  of  time  series  modeling  generally 
provide  a  powerful  method  for  assessing  mean  level 
estimation  precision  and  for  formulating  sampling 
strategies.  Appendix  A  fully  describes  the  use  of  the 
Dynamic  Data  System  (DDS)  for  doing  analyses  based 
on  fitting  Autoregressive-Moving  Average  (ARMA) 
time  series  models  to  the  daily  average  noise  level  data. 
The  analysis  in  Appendix  A  uses  about  1  year's  daily 
CNEL  data  gathered  at  several  sites  around  San  Diego's 
Lindbergh  Field  and  Miramar  Naval  Air  Station  and 
shows  a  high  degree  of  positive  autocorrelation  in 
CNEL  values  from  day  to  day.  This  day-to-day  positive 
autocorrelation  is  to  be  expected  because  of  prevailing 
winds,  slowly  varying  weather  fronts,  and  the  relative¬ 
ly  constant  set  of  daily  operations  and  fleet  mix  at 
commercial  airports.  The  autocorrelated  nature  of  the 
data,  particularly  the  degree  of  positive  autocorrelation 
among  neighboring  observations,  increases  the  amount 
of  consecutive  sampling  required  to  estimate  the  long¬ 
term  mean  level  with  a  given  level  of  precision  over 
sampling  where  independence  is  assumed. 

In  Appendix  B,  the  DDS  method  is  used  to  model 
about  8  months  of  daily  CNEL  data  gathered  at  12 
sites  near  Los  Angeles  International  Airport.  The 
results  given  in  Appendices  A  and  B  form  a  set  of 
guidelines  for  sampling  strategies  near  civilian  airports. 
However,  these  results  use  only  CNEL  data  and  are 
only  for  west  coast  airports. 

Appendix  C  extends  the  analysis  to  east  coast 
airports  and  to  the  use  of  the  DNL  noise  descriptor. 
Specifically,  about  13  months  of  daily  DNL  data 
were  obtained  for  15  sites  near  Boston's  Logan  Air¬ 
port,  and  about  9  months  of  daily  DNL  data  were 
obtained  for  nine  sites  near  the  Washington,  DC, 
Dulles  Airport  and  14  sites  near  the  Washington, 
DC,  National  Airport.  The  DDS  method  was  again 
used  to  model  these  daily  DNL  data.  Monte  Carlo 
simulations  were  performed  to  verify  the  sampling 
requirements  obtained  from  the  DDS  data  modeling 
and  to  study  alternatives  to  consecutive  sampling. 
The  results  of  these  analyses,  along  with  the  results 
given  in  the  other  two  appendices,  are  used  to  form  a 
set  of  guidelines  for  sampling  strategies  in  the  vicinity 
of  civilian  airports. 

The  DDS  method  can  be  used  to  develop  para¬ 
metric  stochastic  time  series  models  of  the  ARMA 
class.  Daily  sound  exposure  (Eq  1),  when  viewed  as  a 

time  series  of  values  X(,  X2 . X^,  has  been 

shown  to  be  well  characterized  by  such  models  (Ap¬ 
pendices  A  and  B).  This  makes  it  possible  to  determine 
the  precision  associated  with  an  estimate  of  the  yearly 


mean  sound  exposure  level  when  the  observed  daily 
values  X,  are  autocorrelated. 

The  general  ARMA  model  is  given  by 

Xt  ~  ^\  Xt  I  +02  X(  2  +  ■  •  •  V*n  Xt  n 
+  a,  ®l  at  i  ®2  at  2  •  •  •  ~  8n  at  m 

[Eql] 

where: 


X,  is  the  noise  level  (daily  average)  for  day  t. 

a,  is  the  random  disturbance  for  day  t. 

01 , ....  0n  are  autoregressive  parameters. 

§l . 0m  are  moving  average  parameters. 

Given  the  time  series  Xt,  the  appropriate  order  of  the 
ARMA  model  (the  proper  values  for  n  and  m)  may  be 
determined  and  the  parameters  of  the  model  may  be 
estimated  by  the  method  of  least  squares.6 

Most  of  the  tilted  ARMA  models  obtained  for  the 
data  analyzed  in  this  report  are  of  relatively  low  order. 
For  example: 

AR  (1):  X,  =  0,  X,_|  +  at 

AR  (2):  X,  =  0,  Xt_j  +  <0 2  xt-2  +  at 


ARMA (l,l):Xj-0j  X,_j  at  —  0]  at _ j 


The  ARMA  models  fit  to  the  daily  sound  exposure 
time  series  Xt  may  be  used  to  estimate  the  precision 
associated  with  the  sample  mean  X.  It  can  be  shown 
(Appendix  B)  that  the  variance  estimate  of  the  sample 


mean  is  given  by 


A 

Variance  (X)  = 


lEq  2] 


*S.  M.  P-^dit  and  S.  M.  Wu,  Time  Series  and  System 
Analysis:  Mounting  Applications  (John  Wiley  and  Sons,  1982). 


where  is  the  residual  mean  square  for  the  fitted 
ARMA  model.  Given  the  above  variance  estimate, 
100  (1  a)  percent  confidence  intervals  of  the  form 

x  -  t  N  (n  +  in),  I  a/2  I  Variance  ( X )  |  1  (hq3| 

may  be  obtained  for  the  true  yearly  mean  sound 
exposure  level,  thereby  providing  an  estimate  of  the 
precision  associated  with  the  sample  mean  X. 

Discussion 

In  summing  the  modeling  results  across  the  Boston 
and  Washington  Airports  (Appendix  C),  it  becomes 
apparent  that  a  number  of  the  sites  exhibit  nonsta¬ 
tionary  behavior;  i.e.,  the  mean  level  changes  over  the 
year.  As  a  result,  long-term  consecutive  sampling 
requirements  are  very  large,  often  constituting  more 
than  one-third  of  a  year.  This  result  is  much  different 
than  the  sampling  requirements  analysis  for  the  west 
coast  airports  which  exhibit,  in  general,  a  stationary 
stochastic  structure  over  an  entire  year’s  data.  In 
attempting  to  delineate  differences  in  the  character¬ 
istics  of  the  east  and  west  coast  airports  and  the 
sampled  data  obtained,  two  observations  can  be 
noted:  (I)  the  west  coast  airports  arc  one-direction, 
onc-runway  airports,  and  (2)  the  monitoring  sites  for 
the  west  coast  airports  are  generally  closer  to  the 
runways.  The  results  relative  lo  the  east  coast  air¬ 
ports  suggest  that  any  analysis  should  be  confined  to 
data  covering  substantially  a  period  of  I  year. 

Figure  5  summarizes  the  results  for  all  of  the  air¬ 
ports  modeled  (airport  sites  only),  including  Dulles, 
National,  and  Logan  from  Appendix  C  and  Los  Angeles 
(including  one  site  from  Lindberg  Field)  from  Ap¬ 
pendices  A  and  B.  This  figure  graphically  represents 
the  similarities  and  differences  among  the  airports  in 
terms  of  their  sampling  characteristics.  The  results 
generally  show  that  the  west  coast  airports  (typically 
one-direction  because  of  prevailing  winds  off  the 
ocean)  tend  to  have  lower  coefficients  of  variation  and 
comparable  autocorrelation  factors  relative  to  the 
multirunway  or  multidirection  (variable  wind)  or  both 
characters  of  east  coast  airports.  These  results  produce 
overall  sampling  requirements  for  the  west  coast 
airports  which  are  generally  lower  than  those  for 
the  east  coast  airports. 

Because  of  the  presence  of  nonstationary  trends 
and  large  sampling  requirements  for  some  of  the  data, 
Monte  Carlo  sampling  experiments  were  performed 
with  the  Los  Angeles.  Boston,  and  Washington  data. 
Through  such  simulations,  generalized  sampling 
strategies,  including  alternatives  to  consecutive 
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sampling,  may  he  examined.  Such  alternate  strategics 
may  require  fewer  total  samples  than  consecutive 
sampling  and  provide  a  means  to  accommodate  trends. 

Sampling  experiments  were  performed  on  those 
sites  with  1  year  of  reasonably  continuous  data.  In 
the  first  set  of  experiments,  the  total  number  of 
samples  taken  to  estimate  the  mean  noise  level  for  each 
strategy  was  28  days;  in  the  second  experiment,  the 
total  number  of  days  was  56. 

Figures  6a  through  6h  show  the  Monte  Carlo 
simulation  for  Los  Angeles,  Boston,  National,  and 
Dulles  airports,  respectively.  For  Los  Angeles,  periodic- 
sampling  indicates  a  slight  but  not  marked  improve¬ 
ment  in  predictive  precision  over  consecutive  sampling. 
With  the  exception  of  two  sites,  the  Los  Angeles  results 
show  that  a  ±  SO  percent  precision  can  be  attained  with 
28  samples,  regardless  of  the  sampling  strategy  chosen. 
For  the  Boston  and  Washington  airports,  significant 
improvements  in  the  predictive  precision  can  be 
achieved  by  periodic  sampling,  e.g.,  1  week  from  each 
quarter  over  the  year.  This  is  particularly  true  for  those 


sites  which  exhibited  nonstationarv  behavior.  To 
guarantee  a  J  60  percent  precision  level  for  these 
airports,  it  is  required  to  sample  1  week  from  each 
quarter  over  the  entire  year. 

In  considering  the  results  of  the  simulation  experi¬ 
ments  involving  requirements  of  56  days  of  sampling, 
it  is  noted  that  for  Los  Angeles.  ±  35  percent  precision 
is  attainable  for  all  sites  regardless  of  the  sampling 
strategy  except  for  Sites  12  and  W3.  For  these  two 
sites  ±  35  percent  precision  is  obtainable  by  sampling 
for  1  week  out  of  each  eighth  of  the  year.  For  the 
Boston  and  Washington  airports,  ±  40  percent  precision 
can  be  achieved  for  all  sites  if  eight  1-week  samples 
are  taken,  one  from  each  eighth  of  the  year. 

For  Los  Angeles,  the  DDS  modeling  consecutive 
sampling  requirements  and  those  obtained  from  the 
Monte  Carlo  simulations  are  generally  about  the  same 
(see  Table  3).  These  sites  exhibit  a  stationary  stochastic 
structure  for  the  entire  year  and  the  simulations  verify 
the  DDS  modeling  results.  The  same  comparison  holds 
for  Boston  sites  1,  5,  8;  Dulles  site  8;  and  National  sites 
1 5  and  20.  These  sites  exhibit  stationary  behavior. 


Table  3 

Comparison  of  DDS  and  Monte  Carlo  Simulation  Results 


Site 

DDS  Modeling  Results 
Consecutive  Samples 

For  P  =  ±50% 

Monte  Carlo  Simulation 
Results  %  P  for 

28  Consecutive  Samples 

Monte  Carlo  Simulation 
Results  %  P  for 

56  Consecutive  Samples 

LAX. 

A1 

6 

26.0 

17.0 

L.A.X. 

A2 

19 

32.0 

22.0 

LAX. 

HI 

4 

24.0 

15.0 

L.A.X. 

12 

8 

32.0 

24.0 

L.A.X. 

12 

30 

77.0 

47.0 

L.A.X. 

L2 

11 

35.0 

18.0 

L.A.X. 

W2 

16 

36.0 

25.0 

L.A.X. 

W3 

60 

81.0 

52.0 

L.A.X. 

W4 

11 

40.0 

28.0 

BOSTON 

1 

11 

40.0 

25.0 

BOSTON 

3 

60 

71.0 

78.0 

BOSTON 

4 

128 

121.0 

116  0 

BOSTON 

5 

35 

65.0 

49.0 

BOSTON 

6 

79 

93.0 

42.0 

BOSTON 

8 

31 

51.0 

39.0 

DULLtS 

1 

90 

57.0 

46.0 

DULLTS 

4 

Nonstationary 

74.0 

47.0 

DULLtS 

6 

158 

93.0 

56.0 

DULLtS 

7 

34 

47.0 

24.0 

Dl'LLIS 

8 

10 

32.0 

20.0 

Dl'LLt’S 

10 

62 

55.0 

39.0 

NATIONAL 

13 

48 

72.0 

46.0 

NATIONAL 

14 

97 

55.0 

52.0 

NATIONAL 

15 

15 

47.0 

35.0 

NATIONAL 

18 

14 

96.0 

900 

NATIONAL 

20 

10 

26.0 

20.0 

NATIONAL 

21 

67 

74.0 

5  7.0 

NATIONAL 

22 

159 

80.0 

66.0 

19 


For  tire  sites  exhibiting  nonslutionury  behavior  at 
Boston,  Dulles,  and  National,  the  comparison  of  the 
DDS  modeling  and  simulation  resulls  are  not  always 
consistent.  In  particular,  the  DDS  models  may  over¬ 
quote  the  consecutive  sampling  requirements  needed 
for  a  particular  level  of  precision  when  compared  with 
the  simulation  results.  This  is  particularly  true  for  site 
6  at  Boston,  site  6  at  Dulles,  and  sites  14  and  22  at 
National. 

Summary 

The  results  generally  show  that  the  west  coast 
airports  tend  to  have  lower  coefficients  of  variation 
and  comparable  autocorrelation  factors  relative  to  the 
multirunway  and/or  multidirection  (variable  wind) 
east  coast  airports.  These  results  produce  overall 
consecutive  sampling  requirements  for  the  west  coast 
airports  which  are  generally  lower  than  those  for  the 
east  coast  airports. 

A  precision  (95  percent  confidence)  of  ±2  to  ±3 
dB  is  generally  achieved  with  4  to  8  weeks  of  mon¬ 
itoring;  1  week  from  each  eighth  or  I  week  from  each 
quarter  of  the  year,  respectively. 
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The  airport  results  (Chapter  3)  show  that  a  preci¬ 
sion  (95  percent  confidence)  of  ±2  to  ±3  dB  is  general¬ 
ly  achieved  with  4  to  8  weeks  of  monitoring;  1  week 
from  each  eighth  or  1  week  from  each  quarter  of  the 
year,  respectively. 

The  point  of  closest  approach  of  airplanes  to  any 
noise  monitor  is  never  more  than  a  few  thousand  feet. 
The  operations  at  civil  aviation  airports  are  quite 
regular  throughout  the  year;  that  is,  the  airport  sched¬ 
ule  remains  relatively  constant.  In  contrast,  training 
operations  at  an  Army  base  are  much  less  regular  and 
the  “point  of  closest  approach"  from  firing  points  to 
monitors  may  be  at  distances  of  5  miles  or  more. 

Since  precision  is  proportional  to  the  square  of  the 
number  of  samples,  an  increase  in  variability  of  only 


a  factor  ol  2  in  the  Army  case  over  the  public  airport 
case  implies  sampling  strategies  lequirmg  1 6-04  weeks 
out  of  the  year  In  effect,  due  to  the  variability  in 
dav  to-day  operations  coupled  with  the  variability 
of  sound  propagation  over  long  distances,  the  only 
way  to  estimate  the  yearly  CDNL  in  the  vicinity  of  an 
Army  base  with  any  degree  of  precision  is  to  measure 
the  CDNL  for  the  entire  ,  ear. 


5  CONCLUSIONS 


1.  The  airport  results  (Chapter  3)  show  that  a 
precision  (95  percent  confidence)  of  ±2  to  ±3  dB  is 
generally  achieved  with  4  to  8  weeks  of  monitoring; 
1  week  from  each  eighth  or  1  week  from  each  quarter 
of  the  year,  respectively. 

2.  The  results  in  Chapter  2  show  that  monitored 
levels  may  differ  substantially  from  computer  simula¬ 
tion  predictions.  In  view  of  the  community  response 
data  which  correlate  well  with  computer  simulation 
and  which  indicate  the  presence  of  substantial  impul¬ 
sive  noise,  it  can  only  be  concluded  that  either  the 
current  monitoring  techniques  are  inadequate  to 
measure  the  true  impulsive  noise  and  that  the  current 
results  are  biased  to  the  low  side  or  that  the  computer 
predictions  are  high  and  communities  respond  ad¬ 
versely  to  much  lower  levels  of  impulsive  noise  than 
is  commonly  believed.  The  former  seems  to  be  the 
more  reasonable  conclusion. 

Because  the  extrapolation  of  the  airport  data 
(Chapter  4)  indicates  impulsive  noise  monitoring 
must  be  continuous  to  properly  measure  the  CDNL 
and  because  survey  data  show  that  even  continuous 
monitoring  will  generally  be  biased  to  the  low  side 
(except  near  base  boundaries),  the  general  recom¬ 
mendation  is  that  monitoring  not  be  performed  until 
monitoring  techniques  can  be  improved  and  a  good 
correlation  achieved  between  monitoring  and  atti- 
tudinal  survey  results.  Until  then,  it  is  recommended 
that  reliance  be  placed  only  on  computer  simulation 
since  these  results  correlate  better  with  attitudinal 
surveys. 


« 


REFERENCES 


Community  Noise  Exposure  Resulting  From  Aircraft 
Operations:  Computer  Program  Description, 
AMR1/TR-73-106  (Department  of  the  Air  Force, 
November  1974). 

DeVor,  Richard  E.,  et  al.,  "Development  of  Temporal 
Sampling  Strategies  for  Monitoring  Nois e,"  Journal 
of  the  Acoustical  Society  of  America,  Volume  66, 
No.  3  (September  I97‘)),  pp  763-771 . 

Guidelines  for  Considering  Noise  in  l.and  Use  Planning 
and  Control  (Federal  Interagency  Committee  on 
Urban  Noise,  June  1980). 

Integrated  Noise  Model  (INM)  (Department  of  Trans¬ 
portation,  Federal  Aviation  Administration.) 

Pandit,  S.  M.,  and  S.  M.  Wu,  Time  Series  and  System 
Analysis:  Modeling  Applications  (John  Wiley  and 
Sons,  1982). 

Pawlowska,  V.,  and  L.  Little ,  The  Blast  Noise  Prediction 
Program:  User  Reference  Manual,  Technical 
Report  N-75/ADA074050  (U.S.  Army  Construc¬ 
tion  Engineering  Research  Laboratory  [CERL] , 
August  1979). 


Schomer,  Paul  D.,  Community  Reaction  to  Impulse 
Noise:  Initial  Army  Survey,  Technical  Report 
N-100/ADA101674  (CERL,  June  1981). 

Schomer,  Paul  D.,  Community  Response  to  Impulse 
Noise:  A  10-Year  Research  Summary,  Technical 
Report  N-167  (CERL,  November  1983). 

Schomer,  Paul  D.,  and  Richard  E.  DeVor,  “Temporal 
Sampling  Requirements  for  Estimation  of  Long- 
Term  Average  Sound  Levels  in  the  Vicinity  of 
Airports,"  Journal  of  the  Acoustical  Society  of 
America,  Volume  69,  No.  3  (March  1981), 
pp  713-719. 

Schomer,  Paul  D.,  et  al..  “Sampling  Strategies  for 
Monitoring  Noise  in  the  Vicinity  of  Airports," 
Journal  of  the  Acoustical  Society  of  America, 
Volume  73,  No.  6  (June  1983),  pp  2041-2050. 

Schomer,  Paul  D.,  et  al..  The  Statistics  of  Amplitude 
and  Spectrum  of  Blasts  Propagated  in  the  Atmos¬ 
phere,  Technical  Report  N-13,  Volume  1  (ADA 
033475)  and  Volume  2  (ADA033361)  (CERL, 
November  1976). 


22 


APPENDIX  A: 

DEVELOPMENT  OF  TEMPORAL  SAMPLING 
STRATEGIES  FOR  MONITORING  NOISE 


R.  E.  DeVor 

University  of  Illinois  at  Urbana  -Champaign.  Urbana,  Illinois  61801 

P.  D.  Schomer 

U.  S.  Army  Corps  of  Engineers,  Construction  Engineering  Research  Laboratory.  Champaign, 

Illinois  6 1820 

W.  A.  Kline 

University  of  tltinoi'-  at  Urbana  -Champaign.  Urbana.  Illinois  61801 

R.  D.  Neathamer 

U.  S  Army  Corps  of  Engineers,  Construction  Engineering  Research  Laboratory,  Champaign, 

Illinois  61820 

(Received  23  November  1977;  accepted  for  publication  12  Apnl  1979) 

Thu  p-vper  addresses  the  problem  of  the  estimation  of  the  long-term  (yearly)  mean  of  the  Community 
Notsc  Equivalent  l  evel  (CNEL)  or  day/night  average  sound  level  (LDs)  Reo-m  environmental  noise 
standards  hwe  emphasized  the  significance  of  this  problem-  While  U  is  possible  to  continually  monitor  (he 
mine  level,  it  is  not  necessarily  desuaMc  or  practical.  It  is  desirable  to  sample  the  level  over  a  relatively 
short  pervxi  of  time  and  use  this  information  to  draw  reliable  inferences  about  the  long  term  mean  level. 
Examination  of  daily  average  noise  levels  (cither  in  mean  square  pressure  or  in  decibel  units)  shows  that 
while  the  data  may  be  stationary  with  respect  to  mean  level  over  a  several  month  period,  tliey  exhibit  a 
strung  pattern  of  autocorrelation  in  which  positive  correlation  predominates.  As  a  result,  the  sample  sires 
required  to  achieve  a  desired  level  of  precision  in  the  sample  mean  estimate  are  much  larger  than  they 
would  otherwise  be  if  the  data  were  uncorrelated  serially  in  tune  To  assess  the  level  of  autocorrelation  in 
the  data,  auioregtcsvi  vc-movmg  average  (AKMA)  models  are  developed  for  the  noise  dau  via  the 
Dynamic  Data  System  (DDS)  approach  to  lime  senes  analysis  These  models  are  then  used  to  derive 
estimates  of  the  sample  mean  variance  and  therefore  to  establish  sampling  strategies.  For  the  data 
examined,  to  obtain  an  estimate  of  the  mean  level  within  a  5-dB  range  (4.509:'  of  the  mean  in  mean 
square  pressure  units),  sample  sizes  in  the  range  of  20-50  consecutive  daily  averages  would  be  required. 

If  the  daily  averages  were  uncorrelated  in  time,  only  5-15  consecutive  daily  averages  would  be  required 
The  data  used  in  this  study  were  obtain  from  continuous  monitoring  at  a  number  of  sites  in  the  vicinity  of 
a  busy  Naval  Air  Station  Some  data  obtained  from  a  Urge  commercial  airport  were  also  analyzed  and 
found  to  have  even  stronger  positive  autocorrelation,  and  therefore  requiring  even  larger  sample  sizes  for 
mean  value  estimation 

PACS  numbers:  43  50  Sr.  4).S0.Qp,  43  50  Lj 


LIST  OF  SYMBOLS 

ith  autoregressive  parameter 

ith  moving  average  parameter 

x, 

A  daily  CNEL  at  time  f 

G, 

ith  Green’s  function  value 

Var(.V) 

The  variance  of  the  sample  mean,  X 

Random  shock  oreurring  at  time  t 

r0 

The  variance  of  the  ,Y,'s 

< 

The  variance  of  the  random  shocks  n, 

N 

The  number  of  observations 

B 

The  backward  shift  operator 

k 

Time  lac 

*l-o  ft 

Unit  normal  statistic  at  the  or  significance 

n 

frtb  lag  covariance  between  A',  and  X,„ 

level 

INTRODUCTION 

Assessment  of  community  noise  Is  an  important  data 
analysis  problem  currently  receiving  considerable  at¬ 
tention  In  both  the  public  and  private  sectors.1, 7  In  this 
paper,  we  will  examine  the  fundamental  problem  of 
estimating  the  long  term  average  noise  level  based  on 
measurements  of  the  Community  Noise  Equivalent 
Level  (CNEL).  Due  to  the  close  correspondence  be¬ 
tween  CNEL  and  the  day  /night  average  sound  level  LUH, 


the  results  herein  should  be  equally  applicable  to  moni¬ 
toring  Lok.  In  this  regard,  the  appropriate  sample  in¬ 
terval  size  and  the  total  amount  of  information  collected 
for  the  specified  sample  interval  are  the  sampling  para¬ 
meters  of  interest. 

To  rationally  embrace  the  problem  posed  above,  it  is 
necessary  to  reveal  the  probabilistic  nature  of  the  noise 
level  signal  as  it  evolves  serially  in  time.  This  is  so 
since  the  techniques  of  estimation  and  statistical  in- 
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ferencc  brought  to  bear  on  such  a  problem  are  depen¬ 
dent  upon  the  underlying  statistical  characteristics  of 
the  data,  in  particular,  the  serial  correlation  of  the 
data.  The  inherent  time-varying  behavior  of  the  data 
suggests  that  its  stochastic  nature  should  be  fully  re¬ 
vealed.  In  fact,  it  becomes  clear  after  only  a  cursory 
examination  of  the  signal  to  be  sampled  that  time  series 
analysis  of  one  form  or  another  will  be  required  to  ob- 
tabi  a  valid  and  efficient  temporal  sampling  strategy. 

In  recent  years,  time  series  modeling  via  the  class  of 
autoregressive- moving  average  (AHMA)  models  has  re¬ 
ceived  considerable  attention  both  in  modeling  strategy 
development  and  modeling  applications,  Recently,  Pan¬ 
dit  and  Wu”  have  proposed  a  new  strategy  for  time 
series  modeling  referred  to  as  Dynamic  Data  System 
(DDS).  The  DOS  approach  combines  the  modeling  of 
deterministic  trends  including  periodicities  by  a  general¬ 
ized  Laplace  transform  and  the  modeling  of  the  remain¬ 
ing  stochastic  variation  including  stochastic  trends  and 
periodicities  by  the  general  class  of  autoregressive- 
moving  average  models.  Since  the  models  developed 
are  parametric  in  nature,  they  have  the  advantage  of 
being  able  to  describe  the  autocorrclatcd  nature  of  the 
data  in  simple,  closed  forms  which  greatly  facilitates 
the  mean  level  estimation  and  inference  problem. 

In  this  paper,  we  will  examine  the  use  of  ARMA 
models  as  an  approach  to  the  characterization  of  com¬ 
munity  noise  level  data.  This  approach  is  considered 
to  have  significant  potential,  since  it  not  only  contri¬ 
butes  to  the  solution  of  the  sampling  strategy  problem 
for  CNEL  mean  level  estimation  and  inference,  but 
also  can  provide  a  methodology  for  monitoring  and  mo¬ 
deling  noise  level  lor  Ihe  dcicction  ol  significant  shifts 
in  level,  and  ran  serve  as  a  basis  for  the  analysis  of  the 
relationship  between  noise  level  and  oilier  physical 
phenomena.  In  the  formulation  of  the  sampling  strate¬ 
gies  presented  here,  an  analysis  is  conducted  on  CNEL 
data  obtained  from  several  sites  in  the  vicinity  of  NAS 
Miramar,  San  Diego,  CA.  The  paper  is  presented  in 
four  main  sections.  The  first  section  provides  a  dis¬ 
cussion  of  the  problems  which  autocorrclatcd  data  pose 
to  the  mean  level  inference  problem  and  formulates  the 
problem  mathematically  via  the  ARMA  class  of  sto¬ 
chastic  models.  The  second  section  briefly  describes 
the  DDS  modeling  methodology,  and  its  application  to 
the  problem  at  hand.  In  the  third  section,  modeling  re¬ 
sults  arc  presented  for  two  noise  time  series  and  an 
assessment  of  the  precision  of  the  estimation  of  the 
mean  noise  fevcl  is  provided.  The  precision  assess¬ 
ment  developed  via  the  DDS  modeling  method  is  com¬ 
pared  to  results  obtained  under  the  assumption  that  the 
noise  data  arc  uncorrclalcd  in  time.  Finally,  the  fourth 
section  summarizes  the  sampling  strategy  requirements 
and  discusses  the  varying  results  over  the  several  sites 
from  which  noise  data  were  obtained  and  analyzed. 

I.  THE  MEAN  LEVEL  INFERENCE  PROBLEM 

At  the  outset  of  this  study,  it  was  proposed  to  formu¬ 
late  a  strategy  for  sampling  the  nuise  level  signal  at  a 
given  location  and  using  the  data  obtained  to  estimate 
the  veorfv  average  noise  level  with  some  prcspecificd 
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level  of  precision.  While  this  information  may  be  quite 
useful  from  an  environmental  impact  point  of  view,  it 
will  be  both  insightful  and  necessary  to  carefully  ex¬ 
amine  (he  time-varying  nature  of  the  signal  over  the 
enure  year.  This  will  not  only  impact  the  estimation 
and  inference  problem,  but  strongly  influence  the  prac¬ 
tical  interpretation  gnrn  to  an  average  or  mean  level 
estim.  <e. 

A.  Autocorrelated  nature  of  the  data 

If  a  scries  of  observations  A',, A', . X„  is  used  to 

estimate  the  average  yearly  noise  level,  the  specific 
nature  ol  the  autocorrelation  in  the  data  will  impact  the 
precision  of  this  estimate.  The  precision  of  the 
estimate  of  the  true  mean  level  is  given  by 


where  y0=Var(A',),  i  =  timc  lag,  N- sample  size,  and 
y,-Fth  lag  auto  covariance  between  X ,  and  A-,.,. 

If  the  data  are  uncorrolatcd,  then  r,  =  0  for  J fc  ;■  1  and 

Var(A')  =  vW.  (2) 

Therefore,  a  (1  -  ar)100?c  confidence  interval  for  X  will 
be  given  by 

X'i*..<,/*(Var(*)r*.  (3) 

Strictly  speaking,  the  I  distribution  should  be  used  in¬ 
stead  of  the  unit  normal,  i  distribution  to  account  for 
the  uncertainty  in  estimating  the  variance  of  the  sam¬ 
ple  mean.  However,  for  the  sample  sizes  encountered 
in  this  paper,  the  l  distribution  is  closely  approximated 
by  the  /  distribution.  While  the  A',  need  not  be  Gaussian, 
they  should  fluctuate  about  a  fixed  mean  level  with  a 
constant  pattern  of  irregularity,  i.e.,  be  stationary. 

If  the  data  are  autocorrclatcd,  then  Eq.  (1)  may  be 
rewritten  as 

Var(A')-  (y„/.V)C ,  (4) 

where  C  is  a  factor  which  varies  with  and  accounts  for 
the  specific  autocorrelated  nature  of  the  data.  In  gene¬ 
ral,  for  positively  correlated  data  the  autocorrelation 
factor  C  will  be  greater  Ilian  1.0  while  for  negatively 
correlated  data  it  will  be  between  0.0  and  1.0.  This 
phenomenon  can  be  appreciated  intuitively  by  examina¬ 
tion  of  Fig.  1.  It  Is  clear  that  for  positively  correlated 
data  the  excursions  or  runs  above  and  below  the  mean 
produce  sample  averages  with  wider  dispersion  about 
the  true  mean  than  if  the  data  were  random.  Similarly, 
negatively  correlated  data  a^e  chai aclenzcd  by  suc¬ 
cessive  high  and  low  values  which  tend  to  "aierage"  to 
values  quite  closely  clustered  about  the  true  mean. 

These  characteristic  behaviors  may  lie  quantified  by 
the  autocorrelation  factor,  which  for  the  modeling 
technique  used  herein  can  be  shown  to  be  solely  a  func¬ 
tion  of  Ihe  parameters  of  an  autoregressive -moving 
average  model  for  the  data. 

B.  Mathematical  statement  of  the  problem 

Consider  a  finite  set  of  discrete  measurements,  .1 , , 
X?,...,XV  obtained  by  uniformly  sampling  3  con'.-, mi  nis 
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(b)  POSITIVE  AUTOCORRELATION 


(C)  NEGATIVE  AUTOCORRELATION 


FIC.  1.  AutocorrcUted  structures  In  time  scries  and  their 
Impact  on  sample  mean  precision. 


signal  A'(/)  at  equispaccd  intervals  A/.  Let  us  assume 
that: 

(1)  The  sample  interval  size  At  is  a  prcspccitied  c<hu 
stant  and  is  of  a  size  sufficient  to  capture  the  structure 
of  the  continuous  signal. 

(2)  vY,}  constitutes  a  stationary  time  series,  i.e., 
fluctuates  about  a  fixed  mean  with  a  constant  pattern 
of  irregularity  and 

(3)  The  length  of  record  A' At  is  sufficient  to  adequately 
encompass  the  significant  long-term  features  of  the  con¬ 
tinuous  signal. 

The  problem  to  be  addressed  is,  then,  to  determine 
the  precision  with  which  the  sample  mean  A'  based  on  ft 
observations  estimates  the  true  mean  level,  or  alter¬ 
natively  lo  specify  the  length  of  record  (number  of  ob¬ 
servations)  needed  to  estimate  the  true  mean  p  by  X 
within  a  certain  prespecified  interval  at  a  given  level  of 
statistical  significance. 

For  a  stationary  time  series,  the  {>,}  converge  to 
zero  as  k  Increases  so  that  for  large  A’,  Kq.  (1)  re¬ 
duces  to  the  approximation, 


The  problem  of  mean  value  inference  then  reduces  to 
the  estimation  of  {?,}. 

We  will  be  concerned  with  the  use  of  parametric  sto¬ 
chastic  models  of  the  autoregressive-moving  average 
(ARMA)  class  lor  the  estimation  ot  the  {>,}.  This  ap¬ 


proach  provides  botli  a  mathematically  sound  and  prac¬ 
tically  appealing  approach  to  the  problem. 

II.  TIME  SERIES  MODELING  BY  DYNAMIC  DATA 
SYSTEM  (DOS) 

A.  The  DDS  modeling  approach 

Modeling  of  stochastic  phenomena  by  the  general  class 
of  autoregressive-moving  average  (ARMA)  models  has 
found  a  tremendous  growth  in  applications  in  the  areas 
of  forecasting  and  control  in  recent  years.  Several 
unified  strategics  have  been  proposed  to  facilitate  this 
modeling  with  varying  philosophies  on  the  model  build¬ 
ing  procedure,  physical  interpretation  of  models,  and 
the  manner  in  which  both  deterministic  and  stochastic 
trends  are  modeled.  The  Dynamic  Data  System  (DDS) 
methodology  has  particular  appeal  for  several  reasons 
which  go  beyond  the  scope  of  this  paper.  In  particular, 
the  modeling  and  interpretation  of  physical  systems  is 
greatly  enhanced  by  this  approach. 

When  only  stochastic  variation  is  evident  in  the  data 
(no  deterministic  trends  such  as  periodicities)  the 
general  class  of  ARMA  models  given  by  Eq.  (6)  is  em¬ 
ployed  for  modeling, 

X , -  +  <t> 2'Y(.2  +  •  •  •  +  ip„A +  — 

“  *  •  •  ”  ®m°t-  at  i  (6) 

where  X ,  Is  the  observation  in  the  time  series  at  time 
/,  a,  is  normally  and  independently  distributed: 

NlD(0,o’),  A, . <\,:  n  autoregressive  parameters, 

and  ,0„:  m  moving  average  parameters.  In  the 

DDS  modeling  methodology,  m  is  generally  defined  to 
be  h-  1  although  the  final  appropriate  fitted  model  may 
have  m  <(n-  1).  It  can  be  shown5  by  employing  the 
elementary  theory  of  linear  operators  on  Hilbert  space 
that  any  stationary  stochastic  system  can  be  approxi¬ 
mated  by  an  autoregressive  moving  average  model  of 
order  («, «  -  1). 

For  the  use  and  interpretation  of  the  ARMA(n,n- 1) 
model  class  it  is  useful  to  consider  two  important  char¬ 
acterizations  of  the  model:  (i)  Green’s  function,  and 
(ii)  autocovariancc  function.  According  to  Wold’s  de¬ 
composition,  X,  may  be  expressed  as  a  sum  of  ortho¬ 
gonal  vectors  G ,a,_,,  in  an  infinite  dimensional  space, 
i.e., 


%G'a< 


The  weights  G,  may  be  determined  for  any  ARMA(n,  n  -  1) 
model  by  equating  coefficients  ot  like  powers  of  the 
backward  shift  operator  B  in. 


(G0-G,R-  C3B’ - )  = 


(1- 0,8-0^'-  . 
(1  -  <f>,8  -  <fi2B! . 


For  stable  systems,  only  a  relatively  few  number  of 
vectors  need  to  hr  added  to  obtain  A",.  The  weights  G, 
are  referred  to  as  Green's  function  and  arc  exprcssable 
in  terms  of  the  tf>’s  and  0’s  of  the  model.  Physically, 
Green’s  function  may  be  thought  of  describing  the  na¬ 
ture  of  the  dynamic  response  of  the  system  to  a  random 


J  Acquit  Soc  Am,,  Vot  68,  No.  3,  September  1979 


DeVor  er  si..  Temporal  sompimg  strategies  for  monitoring  noise 


25 


disturbance  a,.  If  the  forcing  function  a,  is  removed, 
the  system  response  decays  to  the  mean  level  according 
to  t,.  The  Green's  function  rhararten/ution  of  the 
stochastic  process  is  most  useful  in  deriving  statistical 
properties  of  the  process  as  will  be  seen  later. 

The  autocovariance  function  for  the  general 
ARMA(>i,»-  1)  model  form  may  be  derived  from  Fq. 

(6),  noting  that  the  Hh  lag  autocovarianct  y,  is  given 

by 


y,=  E(X,  (9) 

The  fact  that  the  autocovariances  {y,}  can  be  expressed 
solely  in  terms  of  the  model  parameters,  d>.  $.  and  o’ 
is  found  to  be  most  useful  later  when  an  estimate  of  the 
variance  of  the  sample  mean  is  to  be  obtained  given  an 
appropriate  AftMA  model  for  the  data.  The  Appendix 
provides  the  equations  which  define  the  {>,}  as  functions 
of  0,  $,  and  a\. 

In  the  DDS  methodology,  the  appropriate  model  for  a 
given  set  of  data  is  determined  by  successively  fitting 
models  of  progressively  higher  order  by  the  method  of 
least  squares  until  a  satisfactory  fit  is  obtained.  An¬ 
alysis  of  variance  is  performed  for  each  model,  and  the 
y  test  is  employed  to  determine  when  the  reduction  in 
the  residual  sum  of  squares  from  one  model  to  the  next 
is  statistically  significant.  Initially,  an  AItMA(2, 1) 
model,  i.o., 

A ,  -  £|A +  n(  —  (lb) 

Is  fit  to  the  data.  Models  of  the  form  AftMA(4,3), 
ARMA(C,  5),  . . . ,  ARMA(2h,  2u  -  1)  are  then  successively 
lit.  The  model  is  incremented  by  steps  of  two,  i.o., 
ARMA(2h,2h-  1)  so  that  the  roots  of  the  model  can  be 
either  real  or  complex  at  any  time,  thereby  not  forcing 
a  real  root  to  be  present  in  a  model  for  a  process  which 
does  not  physically  have  that  characteristic.  Modeling 
Is  terminated  when  the  />'  test  fails  to  show  significance 
when  the  next  higher-order  model  is  fit.  Individual 
model  parameters  near  zero  may  be  examined  for  sig¬ 
nificance  by  computing  their  (1  -  a)HKTc  confidence 
Intervals.  Insignificant  parameters  are  dropped  and  the 
remaining  parameters  arc  reestimated.  In  general,  an 
ARMA(»,»i)  model  results. 

B.  Estimation  of  the  {yk}  from  the  ARMA  model 
parameters 

By  employing  the  ARMA  model  class  to  characterize 
a  time  scries  of  noise  data,  the  {>,(  for  those  data  (and 
In  Eq.  (5)|  can  be  estimated  by  functions  of  the  model 
parameters  alone.  In  particular, 


Var(.V)  = 

A 


The  variance  of  the  disturbances  a\  may  be  calculated 
by  recursively  calculating  (he  a’ s  from  the  lilted  model 
and  then  substituting  into  Eq.  (12); 


v 


(12) 


The  approach  to  assessing  the  confidence  associated 
with  Hie  sample  mean  .V  of  the  autocorrelated  sequence 
Xp.Y,,. . . , ,YV  is  then  as  follows; 

(1)  Use  the  DDS  modeling  methodology  to  find  by 
successive  fitting  the  appropriate  ARMA(«,m)  model 
for  the  data.  Since,  in  general,  the  models  arc  non¬ 
linear  in  the  parameters,  an  interative  nonlinear  least 
squares  routine  is  required  to  estimate  the  parameters. 

(2)  Based  on  the  fitted  model  and  estimated  parame¬ 
ters  f ,  6,  and  A’,  estimate  the  variance  of  the  sample 
mean  from  Fqs.  (11)  and  (12). 

(3)  Establish  a  (1  -  a)100Cc  confidence  interval  for  the 
true  mean  p  from  Eq.  (3). 

It  should  be  noted  that  in  all  of  the  modeling  ol  noise 
data  which  follows,  the  data  are  modeled  and  precision 
estimates  in  the  sample  mean  are  determined  in  units 
of  sound  exposure  or  mean-square  pressure.  Confi¬ 
dence  intervals  determined  in  mean-square  pressure 
are  then  transformed  to  sound  exposure  level  (SLL)  in 
decibels  by  the  transformation, 

SEL=  101ogin(mean-squarc  pressure).  (13) 


III.  ANALYSIS  OF  NOISE  DATA 

To  illustrate  the  modeling  technique  employed  and  its 
use  in  the  mean  level  estimate  precision  assessment 
problem,  two  sets  of  noise  data  arc  examined.  One  is 
derived  from  a  site  in  the  vicinity  of  NAS  Miramar, 

San  Diego,  CA  (site  30)  while  the  second  was  obtained 
by  nearby  Lindbcrg  Field  (site  W50),  the  commercial 
San  Diego  Airport.  Both  sets  of  data  were  recorded 
between  January  and  June  of  1970.  Figure  2  shows  a 
map  of  NAS  Miramar  and  the  lorations  of  several  moni¬ 
toring  sites  around  the  airfield.  Figures  3  and  4  show 
the  data  in  units  proportional  lo  mean- square  pressure. 
Each  data  point  is  a  time-weighted  24-hour  average 
noise  level  referred  to  as  the  Community  Noise  Equi¬ 
valent  Level  (CNEL).  CNEI,  values  are  determined 
from  the  equation, 

i  r  gas  zoo 

CNEL=  10  lug  l  10 P^Odt 

X19  200  -E0  400  a 

^(fWt-  I  io/*;(/V/  ,  (i4) 

.1200  J|*20 0 

where 20  micropascal  and  r=8G400  s. 

A.  DDS  modeling 

The  DDS  modeling  methodology  was  applied  to  the 
data  to  obtain  an  adequate  AK.MA(n.»i)  model.  For  the 
NAS  Miramar  data  successive  fitting  and  testing  for 
adequacy  via  the  F  test  revealed  that  an  ARMA(8,7) 
model  is  required  to  describe  the  data.  For  the  Lind- 
berg  Field  data,  an  2\1<MA(2, 1)  model  was  found  to  pro¬ 
vide  an  adequate  representation. 

Table  I  provides  the  fitted  models  and  statistical  para- 
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TABLE  I.  Fitted  models  sot  parameter  estimates  for  the  CNF  I.  mean  square  pressure  data  for  sites  HO  and  WSO. 


Site 

Fitted  AHMAOi.m)  model 

X 

>« 

A 

NAS 

Miramar 

site 

Xt  -F  0 .61 6A<  A  -  0 .330.Y,  ♦  0 . 1  VjX,  4 
♦  0 . 11  S.Y,  4  -  0 . 1 3  2A*,  ^  ♦  0  .OOS.Y,  4 

1.60  *  10* 

2.54  *  10” 

1 .75  *  10” 

lor. 

30 

— O.lSliA'i.1  -  0.01 1A*,  nt  ♦  A 

-0.52Haf^+  0.041*1,  .j  ♦ 

-0.119at^+  0.l37al4  -  0.693^^ 

Undl^crg 

Field 

site 

W50 

Xt*0.CMXt4  +  0.135Y,.J 

♦  0.202a, 

St  .9  x  !(,* 

1166.57  x  10” 

750.21  x  10” 

162 

j 

i 

i 

i 


meter  estimates  (or  both  the  site  30  and  site  W50  data. 
As  can  be  seen  from  the  table,  the  two  noise  level  time 
scries  appear  to  vary  considerably  in  terms  o(  both 
their  average  levels  and  stochastic  structure.  For  site 
30,  the  mean  level  is  equivalent  to  61.4  dD  while  for 
site  W50  it  is  considerably  tiigher,  7C.0  dD.  The  dif¬ 
ferences  in  the  autocorreiated  structure  will  be  more 
fully  revealed  when  an  assessment  in  the  precision  of 
the  sample  mean  X'  is  made. 

8.  Mean  level  precision  assessment 

To  assess  the  precision  of  the  estimate  of  the  sample 
mean,  a  (1  -  a)100%  confidence  interval  for  the  true 
mean  p  may  be  determined.  For  the  site  30  data,  using 
Eq.  (11),  and  a*  from  Table  1, 

VSr(X)  =  4.80x  lO10. 

A  95%  confidence  interval  for  n  is  then  given  by  Eq.  (3), 
1.69  x  10*  t  4.29  x  10’. 

In  decibel  units,  the  mean  estimate  is  62.3  dD  and  the 
95%  confidence  interval  is  bounded  by  G1.0  dU  and  63.3 
dB.  The  interpretation  of  this  interval  is  that  we  are 
95%  confident  that  the  true  mean  for  these  data  is  esti¬ 
mated  within  about  i  25%  in  mean-square  pressure  which 
is  about  minus  1.3,  plus  1.0  dD. 

It  is  interesting  to  compare  the  result  above  to  the 
parallel  result  obtained  if  we  assume  that  the  daily 
average  noise  levels  ordered  in  time  are  independent. 

In  this  case,  using  Eq.  (2), 

V5r(X)*2.419x  1010. 

A  95%  confidence  interval  is  given  by 
1.69  x  10%  0.305  x  10*, 

which  says  that  we  can  estimate  the  mean  within  i  18%. 
in  units  of  mean-square  pressure.  Hence,  by  neglect¬ 
ing  the  effect  of  autocorrelation  in  the  data  we  get  the 
impression  that  fur  a  fixed  sample  of  data  (here,  ap¬ 
proximately  one-third  of  a  year)  we  can  estimate  the 
mean  more  precisely  than  what  we  really  can. 

When  analyzing  the  Lindberg  Field  (site  W50)  data, 
the  effect  of  autocorrelation  is  even  more  pronounced. 
Accounting  tor  the  autocorreiated  nature  ot  the  data, 
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Var<X')  =  4876.37  x  10'°. 

Assuming  independence, 

Var(X)  =  652.0  x  1010. 

The  corresponding  95%  confidence  intervals  are  given 
by 

(i)  Assuming  the  data  are  autocorreiated; 

51.9 xl0‘±  6.98 x  10*. 

(if)  Assuming  the  data  arc  independent; 

51.9  x  10®±  2.60  x  10". 

In  other  words,  while  an  uncorrclatcd  analysis  would 
suggest  a  mean  level  estimate  within  i  5%,  based  on  the 
182  data  points,  the  autocorreiated  analysis  shows  that 
our  precision  is  really  only  about  ±  13.4%. 

For  the  site  W50  (Lindberg  Field)  data  it  is  clear  that 
the  effect  of  autocorrelation  in  the  data  is  more  drama¬ 
tic  in  terms  of  the  precision  assessment  of  the  long  run 
sample  mean.  The  specific  effect  in  terms  of  file  form¬ 
ulation  of  sampling  strategics  will  now  be  examined. 

IV.  DEVELOPMENT  OF  SAMPLING  STRATEGIES 
FOR  ESTIMATION  OF  MEAN  NOISE  LEVELS 

A.  Comparision  of  sample  size  requirements: 

Independence  versus  autocorreiated  analysis 

The  question  of  practical  significance  to  the  acousti¬ 
cian  studying  environmental  noise  problems  is:  "How 
long  must  I  monitor  noise  at  a  particular  site  to  be  able 
to  estimate  the  long  run  mean  level  with  a  given  level  of 
precision?"  In  terms  of  flic  data  we  arc  examining  here 
this  question  asks,  "How  many  days  must  be  monitored  ?" 
For  example,  how  large  of  a  sample  of  daily  averages 
must  we  obtain?  The  answer  to  this  question  may  be 
sought  in  an  effort  to  validate  the  predictions  of  a  noise 
level  model  used  in  the  vicinity  ot  a  military  installation 
or  to  assess  and  detect  possible  noise  level  shifts  over 
a  period  of  time. 

Using  the  data  from  several  sites  nr  und  NAS  Mira¬ 
mar  and  vicinity,  we  will  attempt  to  answer  the  above 
question  and  lend  some  insight  to  the  general  |  roblem 
of  sampling  environmental  noise  and  estimating  mean 
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levels  with  some  prcspocificd  level  ol  precision  re¬ 
quired. 

Returning  to  the  site  30  and  site  W  50  dat.i  previously 
modeled,  the  qurlion  is  asked  "How  many  successive 
daily  average  values  would  have  to  be  obtained  to  esti¬ 
mate  the  mean  level  within  «  50  '  in  units  of  mean-square 
pressure’"  Assuming  a  05  I  confidenc  e  interval  of 
sue  <  0.50.Y  (50  i.  precision  in  the  mean  level  estimate) 
is  required,  we  have  that 

t  0.50.V  *  t  2,..  ,;(VirCV)]‘/J .  (15) 

For  the  site  30  data,  answers  to  this  question  are 
compared  assuming  independence  of  the  data  and  ac¬ 
counting  for  the  autocorrelatcd  structure  in  the  data: 

(i)  Independence  analysis 

Combining  Eqs.  (2)  and  (15),  and  solving  (or  A’,  the 
sample  size,  we  find  that  ,V-  12  satisfies  the  pre¬ 
cision  requirement. 

(it)  Autocorrelatcd  analysis 

In  this  case,  combining  Eqs.  (11)  and  (15),  we 
find  that  ,V«23  samples  (daily  averages)  are 
required. 

When  the  site  W50  data  arc  examined  with  respect  to 
this  precision  requirement,  the  cnnip.irtson  is  even 
more  dramatic.  While  an  independence  assumption  pro¬ 
duces  a  sample  sire  require  mom  ot  one  week  (7  days), 
correctly  accounting  for  the  autnrorrclatcd  nature  of  the 
data  produces  a  sample  size  requirement  of  mtr  set  <•« 
weeks  (50  days).  When  the  autocorrelation  is  strong, 
as  is  the  case  (or  the  Lindhcrg  Field  data,  the  effect  of 
neglecting  autocorrelation  in  the  data  is  qu  severe  in 
terms  of  grossly  underestimating  the  length  nl  time  the 
noise  level  must  be  monitored  to  assess  l lie  mean  level 
precision  within  a  given  prosper  died  range. 

B.  Summary  of  modeling  results  and  sample  sire 
requirements 

In  addition  to  the  two  sites  analyzed  above,  data  were 
obtained  from  two  other  sites  (23,31)  in  the  area  of  NAS 
Miramar.  For  one  of  these  .additional  sites  (site  31)  two 


time  series  of  CN’L L  measurements  were  formed  over 
different  lime  periods  to  examine  the  st.ationarity  and 
homogeneity  of  the  data  over  an  extended  period  of  time. 
Figure  2  shows  a  map  of  the  NAS  Miramar  and  the  loca¬ 
tion  of  the  sites  examined.  CNEL  contours  are  also 
shown  on  tins  figure. 

Figures  5  and  6  show  portions  of  thr  data  analyzed  tor 
sites  23  and  31.  Visual  inspection  of  these  two  noise 
time  series  seems  to  indicate  a  marked  difference  in 
the  time-varying  nature  of  the  data.  The  data  for  site 
23  seem  much  more  random  in  nature  than  those  at  site 
31,  which  appears  to  have  more  of  a  positive  correla¬ 
tion  pattern  (almost  a  weekly  pattern)  with  two  large 
"spikes”  around  the  (ourth  and  sixth  weeks. 

Modeling  by  DDS  showed  that  for  the  site  31  data 
(both  noise  time  scries)  an  ARMA(4,3)  model  appeared 
to  adequately  describe  the  autocorrelated  structure  of 
the  data.  For  site  23,  autocorrelated  structure  of  any 
significance  only  seemed  to  appear  at  lags  7  and  14 
(weekly  periodic-type  structure).  A  model  of  the  fol¬ 
lowing  form  was  fit  to  these  data: 

X,  =  -  0.6C7X,. ,  -  0.315X,.,,  +  a,  +  0.355a,., .  (16) 

When  the  parameters  of  the  above  lilted  model  arc  used 
to  determine  the  variance  estimate  of  the  sample  mean 
an  interesting  result  is  obtained.  Contrary  to  the  re¬ 
sults  from  all  other  sites,  the  use  of  Eq.  (11)  to  esti¬ 
mate  Var(X')  versus  Hie  independence  assumption  [Eq. 
(2)|  show  that  smaller  sample  sizes  are  required  when 
we  recognize  autocorrelation  than  when  we  assume  in¬ 
dependence.  In  summary, 

Var(X)  -  10.4  x  10,:  (autocorrelated) , 

Var(A')  -  24.4  x  10,J  (independence) . 

Tins  suggests  that  tins  data  is  predominated  by  negative 
correlation  which  in  fact  explains  the  somewhat  oscilla¬ 
tory  appearance  of  the  data  in  Fig.  5.  It  is  also  noted 
that  the  general  variance  of  this  data  yy,  is  quite  small 
relative  to  the  data  from  the  other  sites  (0.31  x  10‘-  vs 
2.33  and  2.73  x  1015  lor  sites  30  and  31,  respectively). 
Hence,  it  is  not  surprising  that  the  sample  size  rc- 
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FIG.  5.  NAS  Miramar — data 
from  site  23. 
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FIG.  6.  NAS  Miramar — data  from 
site  31. 
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quirements  for  mean  level  precision  assessment  at  site 
23  are  quite  low.  In  fact,  for  a  i  50«  precision  in  the 
mean  level  estimate  (in  units  of  mean-square  pressure) 
only  one  daily  average  is  required  (two  daily  averages 
if  autocorrelation  is  neglected).  To  obtain  a  ±  23x  pre¬ 
cision,  only  three  daily  readings  are  required  (seven 
daily  readings  if  autocorrelation  is  neglected). 

Table  11  summarizes  the  modeling  and  sampling  stra¬ 
tegy  analysis  for  all  of  the  sites  evaluated.  It  is  inte¬ 
resting  to  observe  the  wide  variation  in  the  autocorre- 
laled  nature  of  the  data  from  these  four  sites  and  its 
concommittant  impact  on  the  sample  size  requirements. 
The  following  observations  are  summarized: 

(1)  Sites  30  and  31  (Figs.  3  and  6)  arc  typified  by  data 
which  have  two  structural  components  visible  in  the 
time  scries:  (i)  irregular  runs  of  average  length  of 
about  7  days,  and  (ii)  a  few  very  sharp  spikes  indicating 
that  high  mean  square  pressure  levels  occur  somewhat 
infrequently  over  an  extended  period  of  time.  Table  11 
indicates  that  the  time  series  modeling  approach  re¬ 
sponded  to  these  data  characteristics  by  conveying  a 
general  positive  correlation  content  which  produced 
sample  size  requirements  in  access  of  those  required 
under  the  assumption  of  no  correlation  in  the  data. 


TABLE  II.  Summary  of  sample  size  requirements. 


Sample  size  requirements  are  about  100%  greater  than 
would  be  called  for  assuming  independence  of  the  data. 

(2)  For  site  W50  (Lindberg  Field),  Fig.  4  shows  the 
same  general  patterns  as  sites  30  and  31  except  that  the 
irregular  runs  (above,  below  the  mean)  seem  longer. 
This  suggests  even  stronger  positive  correlation.  The 
DDK  model  for  this  site  confirms  this  general  appear¬ 
ance  by  producing  very  high  sample  size  requirements 
relative  to  the  assumption  of  independence  (50  days  ver¬ 
sus  7  days). 

(3)  For  site  23  (Fig.  5)  the  data  arc  quite  uniform  in 
variation  level  (no  large  spikes)  and  have  much  more 
of  an  oscillatory  pattern  as  opposed  to  a  pattern  of 
longer,  irregular  runs.  The  DDS  model  conveys  this 
mathematically  by  producing  a  sample  size  requirement 
which  is:  (i)  much  lower  than  for  the  other  sites,  and 
(ii)  such  that  actually  fewer  samples  arc  required,  rel¬ 
ative  to  assuming  independence. 

V.  CONCLUSIONS 

In  this  paper,  a  method  for  developing  strategies  for 
temporal  sampling  of  environmental  noise  has  been  pre¬ 
sented.  Data  were  examined  from  several  sites  in  the 
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arc*  of  NAS  Miramar  and  a  site  from  l.indberg  Field, 
San  Diego.  These  data  are  typified  by  a  strong  auto- 
correlated  structure.  This  autocorrclated  structure 
greatly  influences  the  estimation  of  the  variance  of  the 
sample  mean  and  must  be  accounted  for  in  the  develop¬ 
ment  of  valid  sampling  strategies.  The  Dynamic  Data 
System  approach  was  used  to  develop  stochastic  auto¬ 
regressive-moving  average  models  capable  of  cap¬ 
turing  this  structure  in  a  closed  parametric  form.  The 
parameters  of  these  models  are  used  in  the  estimation 
of  the  sample  mean  variance.  While  the  specific  sampl¬ 
ing  plans  obtained  may  be  peculiar  only  to  the  airports 
studied,  the  modeling  and  inference  techniques  presented 
may  be  applied  to  a  wide  range  of  mean  level  inference 
problems  when  data  arc  autocorrelatcd. 

For  the  data  examined,  comparisons  of  the  sampling 
requirements  assuming  independence  and  accounting  for 
autocorrelation  show  that,  in  general,  more  samples 
are  required  accounting  for  autocorrelation  due  to  the 
predominance  of  positive  autocorrelation  in  the  data. 

By  neglecting  the  effect  of  autocorrelation  in  the  data, 
the  impression  is  given  that  for  a  fixed  sample  of  data, 
one  can  estimate  the  mean  more  precisely  than  what 
one  really  can.  For  example,  when  the  site  W50  data 
are  examined  the  independence  assumption  produces  a 
sample  size  requirement  of  one  week  (7  days),  while 
correctly  accounting  for  the  autocorrelatcd  nature  of 
the  data  produces  a  sample  size  requirement  of  over 
seven  weeks  (50  days). 

Wide  differences  in  the  autocorrelation  structure  and 
the  sampling  requirements  from  site  lo  site  at  the  same 
airjiorf,  NAS  Miramar,  were  noted.  Further  work  is 
being  done  ta  study  the  effects  of  weather,  flight  pat¬ 
terns,  and  monitor  location  on  the  problem  of  mean 
sound  level  assessment. 
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Community  noise  temporal  sampling  requirements  in  general,  and  in  the  vicinity  of  airports  or  other  large 
noise  producers  m  particular,  are  not  well  understood  Frequently,  the  purpose  is  to  sample  and  estimate  the 
true  yearly  Day/Night  Average  Sound  Level  |DNL|  or  Community  Noise  Equivalent  Level  (CNELl  This 
being  the  case,  it  is  important  to  note  that  day-to-day  samples  arc  not  independent,  but,  in  fact,  the  time  senes 
formed  by  these  day-to-day  samples  exhibits  an  autocorreiated  structure.  For  this  reason,  sampling 
requirements  are  4  to  8  times  larger  than  are  calculated  by  assuming  purely  random  day-to-day  data 
Moreover,  the  data  may  exhibit  weekly  and  yearly  deterministic  trends  As  a  result  of  those  factors,  the 
analysis  herein  shows  that  sampling  requirements  sufficient  to  achieve  a  precision  of  *  2  to  -  3  dB  of  the 
true  yearly  CNEL  or  DNL  value  with  a  93%  confidence  level  can  be  summarized  as  1 4  days  of  totally 
random  sampling  throughout  the  year,  or  3  -4  weeks  of  quasi-random  sampling  taken  one  week  at  a  time,  or 
at  least  30  days  of  totally  continuous  sampling 

PACS  numbers  43  30  Lj.  43  30  Sr.  43  30  Nm 


INTRODUCTION 

Recently,  increasing  attention  has  been  given  to  the 
problems  associated  with  high  sound  exposure  levels 
in  the  Immediate  vicinity  of  installations  such  as  civil 
airports  and  military  bases.  From  an  acoustical  point 
of  view,  work  Is  proceeding  on  several  fronts  In  th« 
areas  of  Improved  equipment  design,  better  operations 
planning. and  new  techniques  for  nolae  abatement.  An 
Important  element  of  the  overall  problem  Is  the  mea¬ 
surement  of  nolee  levels  and  the  associated  statistical 
assessment  of  the  precision  of  mean  level  estimates. 
Most  techniques  In  use  today1' 1  for  sampling  community 
noiae  call  for  sampling  over  relatively  short  periods 
of  time.  e. g  from  a  few  minutes  to  perhaps  a  single 

day.  However,  the  time  varying  nature  of  noise  data 
when  viewed  as  a  time  series  (hourly  or  dally  averages) 
suggests  that  short-term  sampling  may  lead  to  serloue 
Inaccuracies  in  the  estimation  of  a  long-term  (yearly) 
average  noise  level.  For  example,  the  24-h  periodic 
pattern  In  hourly  average  sound  level  may  vary  from 
about  40  to  85  dB.  The  Community  Noise  Equivalent 
Level  (CNEL)  or  Day  /Night  Average  Sound  Level  (DNL) 
bath  commonly  vary  from  45  to  80  dB.  These  wide 
ranges  for  sound  level,  together  with  the  fact  that  the 
data  tn  general,  exhibit  high  positive  autocorrelation 
and  high  coefficients  of  variation  suggests  that  small 
and/or  short  sampling  periods  may  provide  both  Im¬ 
precise  and  Inaccurate  mean  value  estimates. 

The  techniques  of  lime  series  modsllng,  tn  gsnsral, 
provide  a  powerful  methodology  tor  an  assessment  of 
mean  level  estimation  precision  snd  the  formulation  of 
sampling  strategies.  In  a  recent  paper,’  the  authors 
have  described  the  use  of  the  Dynamic  Data  System 
(DOS)  for  performing  these  analyses  based  on  the  fitting 
of  Autoregressive -Moving  Average  (ARMA)tlme  series 
models  to  the  daily  average  noise  level  data.  The  pre¬ 
vious  paper  utlUted  approximately  one  year's  dally 


CNEL  data  gathered  at  several  sites  around  San  Diego's 
Lindbergh  Field  and  Miramar  Naval  Air  Station.  Anal¬ 
ysis  of  these  data  ahowed  a  high  degree  of  positive  cor¬ 
relation  in  CNEL  values  from  day  to  day.  This  auto- 
corralated  nature  of  the  data,  particularly  the  degree 
of  positive  correlation  between  neighboring  observa¬ 
tions,  Increases  the  amount  of  consecutive  sampling 
required  to  estimate  the  long-term  mean  level  with  a 
given  level  of  precision  over  sampling  where  indepen¬ 
dence  is  assumed. 

In  this  paper,  the  dynamic  data  system  method  is 
used  to  model  approximately  eight  months  of  daily 
CNEL  data  gathered  at  12  sites  in  the  vicinity  of  the 
Los  Angeles  International  Airport  (LAX).  The  results 
of  this  analysis,  along  with  the  previous  results,  are 
used  to  form  a  set  of  guidelines  for  sampling  strategies 
In  the  vicinity  of  airports.  Of  necessity,  this  paper 
uses  CNEL  data.  However,  as  stated  In  several  ref- 
erences,*’*  the  correspondence  and  correlation  between 
CNEL  and  DNL  are  so  close  as  to  make  the  results 
equally  applicable  to  DNL  data. 

i.  THE  DATA  BASE 

Continuous  daily  monitored  CNEL  values  were  sup¬ 
plied  by  LAX  for  the  period  from  1  May  1977  to 
31  December  1977  for  the  12  monitoring  locations  in¬ 
dicated  in  Fig.  1.  These  daily  CNEL  were  transformed 
to  values  denoted  by  X,  on  a  linear  scale,  proportional 
to  dally  sound  exposure  (SE)  as  shown  in  (1). 

V  JQCSILM0 

Figure  2  illustrates  the  resulting  data  plotted  as  a 
function  of  day  of  the  year  for  these  12  sites  after 
transforming  ihe  data  as  described  in  Eq.  1.  II  musi 
be  noted  (hat  these  V  values  are  used  because  one  is 
interested  in  estimating  the  yearly  mean  DNL  or  CNF1. 
value  which  by  definition  is  estimated  by  the  sample 
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yearly  moan  (arithmetic  average),  V.  Monitored  data 
is  frequently  measured  for  lur  less  lhan  one  year's 
1 1 mo  period,  and  it  is  the  purpose  of  (his  paper  to  ad¬ 
dress  On*  validity  of  techniques  to  estimate  this  yearly 
moan. 

In  order  to  draw  the  most  complete  conclusions  pos  - 
slide,  the  data  developed  in  the  previous  paper  will  be 
utilized  tor  the  analysis  herein  alone  with  the  I, AX 
data.  Figure  3  illustrates  the  location  of  two  monitor¬ 
ing  points  utilized  in  the  vicinity  of  Miramar  Naval  Air 
Station.  Other  sites  at  Miramar  Naval  Air  Station  were 
either  outside  the  noise  contour  area  and  measuring 
other  noise  (e.g.,  site  23)  or  they  had  significant  gaps 
m  the  data  when  ordered  serially.  A  third  set  of  pre¬ 
vious  data  comes  from  a  monitor  located  approximate¬ 
ly  one  mile  west  of  the  main  runway  at  San  Diego's 
Lindbergh  Field. 


Using  the  DDS  method,  ARMA  models  of  various 
orders  were  developed  for  15  of  (he  16  above  lime 
series  of  the  quantity  X.  For  details  of  the  modeling 
methods  and  terminology.  I  he  leader  is  referred  to  tin* 
author’s  previous  paper.  '  One  of  the  LAX  time  senes 
(II,  the  location  of  which  is  shown  in  Fig.  1)  was  found 
to  contain  strong  deterministic  trends  and  thus  was  cut 
modeled  using  the  DDS  stochastic  models.  Using  the 
method  described  in  the  previous  paper,  the  estimated 
parameters  of  the  fitted  ARMA  models  were  employed 
to  determine  the  autocorrelation  factors  lor  each  « •  f 
the  time  senes  and  thereafter  estimate  the  variance  <if 
the  sample  means.  Table  I  summarizes  the  model 
type  and  autocorrelation  lartor  for  each  of  the  time 
series  modeled.  This  Table  also  contains  the  sample 
mean  and  sample  variance  derived  trom  these  time 
senes.  To  aid  in  interpret  at  ion ,  the  coefficient  of 
variation,  which  is  the  ratio  of  the  standard  deviation 
to  the  mean,  is  also  included  in  Table  I. 

Based  on  the  assumption  of  independence  of  the  data 
(no  autocorrelation  on  a  day-to-day  basis),  one  can 
calculate  the  number  of  samples  required  to  estimate 
the  long-term  (yearly)  mean  for  any  desired  level  of 
precision.  In  this  paper,  an  estimation  precision  of 
t  StfA  of  the  mean  with  a  95'.,  confidence  level  has  been 
chosen  as  typical  values  for  illustrative  purposes.  It 
must  be  noted  that  a  plus -minus  50\  band  in  the  esti 
mat) on  of  mean  X  corresponds  to  a  plus  2-  to  minus 
3-dB  band  on  the  estimation  of  yearly  mean  DNL  or 
CNEL.  However,  the  fitting  of  an  ARMA  model  to  a 
time  series  indicates  that  the  series  possesses  a  posi¬ 
tive  autocorrelative  structure.  And  hence,  sample 
sizes  determined  assuming  independence  will  under¬ 
estimate,  in  some  cases,  by  a  considerable  amount, 
the  actual  sample  size  requirements.  When  the  auto¬ 
correlation  factors  are  rightfully  applied,  the  correct 
sample  size  requirements  emerge.  Table  I  also  pro¬ 
vides  the  sample  size  requirements  for  the  estimation 


TABLE  I.  Modeling  results,  summary  statistics,  and  sample  size  requirements  for  each  site. 


■ 

Coefficient 

Sample  size  * 

Sample  size  J 

Mean 

of 

Autocorrelation 

requirement 

requirement 

Site 

Mode! 

X 

Variance 

variation 

factor 

(Independence) 

(autocorrelation) 

ARMLA(2.  1) 

3.89*  101 

2.48*10'* 

ARMA(S.S) 

6  80*  10‘ 

8  38*10" 

AR<1) 

3.88  *10’ 

1.17*10" 

ARMAW.3) 

1.03*10' 

8.70*10" 

AR(1) 

7.03*10' 

8  48*10" 

AR(1) 

2. 87*  10" 

1.03*10" 

ARMA(2, 1) 

2  21  *  10' 

5.14*10" 

Random 

7.44  *  10* 

5.71  *  10" 

ARMA(2,  1) 

8.90*  10T 

6  85*  10** 

AHMAU.l) 

5  03*  10' 

6.17*10" 

ARMA(H,  7) 

1  84*10* 

2.63*10" 

ARM  A  (4, 3) 

1.58*  10* 

2.33*  10" 

AH  MA(2,  1) 

0.19*10' 

1.19*  lO" 

Sample  size  refers  to  the  number  of  consecutive  sample  days  needed  to  predict  the  long-term  mean  level  within  plus  2  to  mmu« 
:<  dH  of  the  true  value  with  95%  confidence  based  on  the  quantity  X  In  Eq.  1.  (Independence)  means  sample  sizes  baaed  on  the 
assumption  of  serially  Independent  data.  (Autocorrelation)  means  sample  sizes  based  on  the  assumption  of  serially  sutocor- 
related  data. 

s  Sites  at  Miramar  Naval  Air  Station. 

‘  Site  at  Lindbergh  Field,  San  Diego,  California. 
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of  the  mean  noise  level  for  both  the  cases  of  assumed 
independence  and  correctly  accounting  for  the  autocor- 
related  structure  in  the  data.  In  the  Table,  the  auto¬ 
correlation  factor  relationship  is  not  precisely  evident 
because  the  sample  numbers  have  always  been  rounded 
up  to  the  next  highest  integer  in  order  to  guarantee  the 
stated  precision.  Table  I  also  provides  Summary  sta¬ 
tistics  and  sample  size  requirements  for  the  two  moni¬ 
tor  sites  at  Miramar  Naval  Air  Station  and  the  one 
site  at  Lindbergh  Field.  These  results  were  previous¬ 
ly  reported  in  Ref.  3  and  are  included  here  to  demon¬ 
strate  the  similarity  in  the  results  across  the  three 
airports. 


In  summary.  Table  I  lists  the  monitor  sites,  the 
model  type,  the  mean,  and  the  standard  deviation  for 
the  original  time  series,  the  coefficient  of  variation, 
the  number  of  independent  samples  required  for  t50'( 
accuracy  (»2  to  -3  dB),  the  autocorrelation  factor  and 
the  true  number  of  samples  required  for  *5(H  accuracy 
when  the  autocorrelated  nature  of  the  time  series  is 
taken  into  account. 


Operational  data  were  also  supplied  by  the  Los 
Angeles  Department  of  Airports  for  landings  per  day  by 
runways  on  25R,  25L;  24R  and  24L;  and  by  runway  pairs 
for  6L  and  6R;  and  for  7L  and  7R.  Takeoff  data  were 
supplied  per  day  by  runway  pair  for  25L  and  25R,  24 L 
and  24R,  6L  and  6R,  and  for  7L  and  7R.  Operations  at 
LAX,  Miramar,  and  Lindbergh  are  typically  westward 
due  to  prevailing  winds  off  the  ocean.  Occasionally 
winds  arc  such  that  the  normal  direction  of  operations 
must  reverse  and  takeoffs  and  landings  are  to  the  east. 


pairs  were  developed  and  are  listed  in  Table  11  In  each 
case,  correlation  coefficients  (i.e.,  the  zeroth  lag 
cross -correlations)  between  the  X  lime  series  and  the 
corresponding  runway  operations  time  series  in  terms 
of  total  daily  approaches,  departures,  and  the  sum  of 
approaches  and  departures  were  determined.  As  can 
be  seen  from  the  descriptions  in  the  Table,  data  pairs 
have  been  selected  to  emphasize  the  predominant  type 
of  operation  likely  to  be  encountered  at  any  given  mon¬ 
itoring  location.  For  example,  because  of  the  west¬ 
bound  nature  to  the  traffic  flow,  site  L2  should  pre¬ 
dominantly  measure  landings  on  25L  and  location  Y.2 
should  measure  both  takeoffs  and  landings  in  either  di¬ 
rection  on  the  south  complex  (where  “complex"  defines 
either  the  northern  or  southern  pair  of  runways  at 
LAX).  Table  fl  provides  only  those  estimated  cross - 
correlations  which  were  deemed  significant  in  magni¬ 
tude.  While  many  site  noise -operations  pairs  were 
correlated,  those  not  found  in  the  Table  resulted  in 
small  (effectively  zero)  correlation  coefficients.  Based 
on  the  amount  of  data  available,  correlation  coefficients 
in  excess  of  approximately  plus  or  minus  0.15  would 
be  considered  statistically  significant  at  the  0.05  level 
of  significance.  However,  for  purposes  of  explaining 
the  data  and  analysis  herein,  only  the  major  significant 
correlations  are  presented  and  discussed.  As  explained 
later  in  this  paper,  0.3  is  chosen  as  the  value  to  define 
major  significant  correlations. 


II.  DISCUSSION  OF  RESULTS 


To  lest  the  relation  of  the  monitored  data  in  the  vi¬ 
cinity  of  LAX  with  operations,  various  correlation 


Examination  of  the  data  in  Table  I  shows  a  wide  range 
of  sampling  requirements  from  site  to  site  depending  on 
relative  location  and  proximity  to  the  runways.  In 
many  cases,  the  sample  size  requirements  are  quite 


TAIll.E  II.  Zeroth  lug  cross- correlation  between  noise  level  recorded  at  a  site  and  operations  In 
Ihe  vicinity  of  the  Kite. 


Operations  strongly  correlated  with 


Correlation  coefficient 


M 

None 

•  .  • 

A2 

EB/APP/7  •  WB/DEP/25* 

0.481 

El 

WB/APP/25 

0.308 

WB/DEP/25 

0.312 

E2 

WB/APP/25 

0.408 

WB/DEP/25 

0.323 

11 

WB/APP/24R1’ 

0.731 

WB/APP/24 

0.639 

WB/APP/24  *  EB/DEP/6 

0.611 

12 

None 

LI 

WB/APP/25 

0.436 

WB/APP/25*  EB/DEP/7 

0.406 

12 

WB/APP/25L 

0.716 

Wl 

None 

... 

W2 

Non* 

W3 

WB/APP/24R  b 

WB/APP/24 

0.499 

WB/APP/24  •  EB/DEP/6 

0.439 

W4 

WB/APP/24Rh 

0.419 

WB/APP/24 

0.432 

WB/APP/24  *  EB/DEP/6 

0.473 

*  EM/AI*l'/7  -  W M/DE  P/20  denotes  the  sum  of  eaathound  approaches  on  the  7  complex  and  wcsi^ 
I* mud  departure*  on  the  20  complex. 

;•  111  dciiotfN  <he  24  complex,  right  runway 
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large.  The  large  sample  numbers  may  be  due  to  the 
presence  of  strong  positive  autocorrelation,  targe 
overall  noise  series  variability  from  day  to  day,  or 
both.  One  could  hope  to  observe  that  the  number  of 
autocorrelated  sample  days  required  grows  smaller  as 
the  sample  site  approaches  the  airport.  Unfortunately, 
this  is  not  evident  since  stations  A2,  W2,  and  W3  each 
require  substantial  numbers  of  days.  However,  care¬ 
ful  examination  of  these  data,  the  operational  data,  and 
weather  conditions  indicates  possible  explanations  for 
the  greater  variability  found  in  two  of  these  three  sta¬ 
tions. 

Site  W3  has  an  autocorrelation  factor  of  4.5,  which  is 
typical  of  the  values  found  lor  all  the  monitoring  sites 
in  close  proximity  to  the  runways  (except  for  A2  and 
W2  which  are  the  other  two  sites  near  the  runways 
which  have  relatively  large  sampling  requirements) 

Site  W3,  however,  exhibits  a  much  higher  coefficient  of 
variation  than  do  the  other  sites  in  the  vicinity  of  the 
airport  (except  for  site  W2  which  is  again  one  of  the 
three  sites  under  discussion  and  site  12  which  will  be 
discussed  later).  Examination  of  the  number  of  land¬ 
ings  per  day  on  runway  24R,  the  operations  which  most 
influence  the  noise  received  at  Site  W3,  indicates  a 
high  degree  of  variation  from  day  to  day.  The  corre¬ 
lation  coefficient  between  the  landings  on  runway  24R 
and  the  noise  measured  at  Site  W3  is  0.673.  Although 
not  shown  herein,  the  number  of  landings  per  day  on 
24R  are  much  more  variable  than  are  the  numbers  of 
landings  per  day  on  the  other  three  runways  (24L,  24R, 
and  25L).  Thus  the  high  coefficient  of  variation  in  the 
X  data  from  Site  W3,  in  reality,  likely  reflects  the  high 
variation  in  the  operations  data. 

Site  A2  exhibits  a  coefficient  of  variation  which  is 
much  in  line  with  the  other  stations  near  to  the  run¬ 
ways  (except  W2  and  W3).  At  this  site,  however,  the 
autocorrelation  (actor  is  8.2— a  value  which  is  much 
higher  than  the  value  for  most  of  the  other  sites.  Ex¬ 
amination  of  this  time  series  shows  that  the  noise  level 
generally  rises  during  the  warmer  summer  months. 
Typically  at  Los  Angeles,  there  are  12  days  in  which 
the  high  temperature  exceeds  90°F.*  This  site  mea¬ 
sures  predominantly  takeoff  noise  and  Is  some  4000  ft 
further  from  start  of  roll  than  Is  site  Al.  Thus  this 
site  is  the  only  site  which  will  be  Influenced  by  the 
average  temperature  since  temperature  strongly  affects 
the  efficiency  of  the  turbojet  engines— requiring  longer 
takeoff  rolls  and  lower  altitudes  over  site  A2  during 
warm  weather.  In  fact,  this  is  the  exact  trend  strongly 
evident  in  the  data;  l.e.,  the  average  sound  level  goeB  up 
during  the  warm  weather  months.  Referring  to  Table 
II,  this  conclusion  Is  further  supported  by  the  fact  that 
the  cross-correlation  between  noise  level  monitored  at 
site  A2  and  the  combination  of  eastbound  approaches 
and  westbound  departures  on  the  south  complex  is  quite 
high,  a  value  of  0.481. 

Site  A2  exhibits  two  predominant  peaks  in  the  data  in 
addition  to  the  general  trend  discussed  above.  One 
peak  appears  on  the  10th  of  May  (day  10)  and  the  other 
on  the  16th  of  December  (day  230).  Many  of  the  sta¬ 
tions  exhibit  peaks  in  the  nild-to-end  of  December  lime 


period.  This  time  period  was  a  period  of  heavy  cloud 
cover,  variable  winds,  and  rain.7  No  explanation  can 
be  found  in  the  weather  or  the  operations  for  the  peak 
exhibited  on  May  10.  Similarly,  examination  of  the 
data  at  site  W2,  the  other  “problem  sue,’  exhibits  a 
very  strong  peak  on  the  16th  of  duly  (day  77).  Again, 
no  weather-related  or  operations-related  explanation 
can  be  found.  The  two  series  of  quantity  .V  values  de¬ 
rived  from  the  monitored  CNEI.  values  at  these  two 
sites  were  remodeled  with  the  data  for  July  16  deleted 
from  the  series  at  site  W2  and  the  data  for  May  10  de¬ 
leted  at  site  A2,  respectively. 

At  site  W2,  removal  of  the  spike  on  July  16  did  not 
change  the  model.  That  is,  the  site  noise  series  re¬ 
mained  totally  random  and  removal  of  the  spike  only 
served  to  decrease  the  variance,  changing  the  sampling 
requirement  from  20  to  16  days.  Similarly,  al  site  A2, 
removal  of  the  spike  on  May  10  did  little  fo  change  its 
results.  However,  removing  both  the  spike  on  May  10 
and  the  spike  on  December  16  substantially  altered  the 
characteristics  of  the  time  senes  hut  did  little  to 
change  the  ultimate  sampling  requirements.  Removal 
of  these  two  spikes  revealed  a  model  with  a  strong  de¬ 
terministic  trend.  The  time  series  strongly  demon¬ 
strated  a  7 -day  weekly  cycle  to  the  data  in  addition  to 
the  yearly  cycle  discussed  above.  Because  of  this 
strong  weekly  cycle,  the  autocorrelation  factor  actually 
rose  from  its  already  high  value  of  8  to  almost  25. 
Similar  results  were  found  in  the  author’s  previous 
paper.1  However,  removal  of  the  spike  at  December 
18  may  not  be  justified  because  its  presence  can  be 
explained  by  weather-related  factors. 

Examination  of  the  cross-correlation  data  between 
operations  and  measured  noise  levels  reveals  some 
correlations  to  be  as  expected  and  others  to  present 
some  significant  departures  from  expectations.  The 
correlations  between  noise  levels  in  the  vicinity  of  the 
east  end  of  the  north  complex  (24),  and  the  correspond¬ 
ing  north  complex  operations  exhibit  the  most  regulari¬ 
ty  and  are  most  as  expected.  That  is,  the  landings  on 
24R  are  highly  correlated  with  the  data  measured  at 
W3,  the  landings  on  24L  are  well  correlated  with  data 
measured  at  site  W4.  Also  the  data  measured  at  these 
two  sites  correlate  well  with  the  total  operations  over 
the  east  end  of  the  north  complex.  The  correlations 
with  the  data  measured  at  site  II,  indicate  that  landings 
on  24R  correlate  with  these  measurements  but  not  the 
landings  on  24L.  This  also  is  a  reasonable  result.  It 
Is  noted  that  the  cross-correlations  between  operations 
on  the  24  complex  and  site  W2  (to  the  north-side  of  the 
complex)  are  generally  small  (effectively  zero).  This 
site  Is  typified  statistically  as  a  random  noise  series 
with  high  coefficient  of  variability  (1.05). 

Unfortunately,  landings  on  the  25  complex  do  not  ex¬ 
hibit  the  same  regularity  and  expected  results  as  de¬ 
scribed  above  (or  the  24  complex.  Landings  on  25L 
are  well  correlated  with  the  noise  measured  at  site  L2, 
but  landings  on  25R  were  not  found  to  be  correlated 
with  the  data  measured  at  site  1.1.  However,  (lie  total 
of  all  operations  or  Ihe  tolal  number  of  landings  to  the 
east  of  the  south  complex  are  botti  well  correlated  wit li 
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the  data  measured  at  site  LI.  This  would  seem  to  in¬ 
dicate  that  site  LI  is  such  that  it  is  located  more  nearly 
acoustically  midway  between  the  operations  on  25 L  and 
25H,  rather  than  as  indicated  in  Fin.  1.  On  the  other 
hand,  the  lower  correlation  for  site  L2  with  overall 
operations  and  the  generally  higher  correlation  with  the 
specific  landing  operations  on  25L  indicate  that  it  is 
more  nearly  m  line  with  25 L.  Station  12  exhibits  no 
correlation  with  any  of  the  operations  either  taken 
singly  or  in  combination  and  like  site  W2,  its  noise 
series  exhibits  weak  autocorrelation  and  u  high  coef¬ 
ficient  of  variation. 

It  must  be  noted  that  correlations  developed  between 
numbers  of  operations  and  measured  data  cannot  lie 
expected  to  be  extremely  high  because  the  measured 
data  are  being  correlated  with  operations  which  may 
represent  a  variety  of  noise  levels.  For  example,  cor¬ 
relations  have  been  calculated  with  landings  alone  when 
in  reality  on  certain  days  the  landings  may  be  very  low 
and  the  takeoffs  very  high  with  the  resulting  noise  levels 
also  high.  On  the  other  hand,  correlations  have  been 
developed  with  total  operations  where  there  is  no  guar¬ 
antee  that  these  operations  do  not  produce  systematical¬ 
ly  differing  CNEL  on  differing  days  that  is  not  re¬ 
flected  merely  in  the  total  number  of  operations.  As  a 
result,  correlation  coefficients  in  excess  of  approxi¬ 
mately  0.3  are  considered  significant  at  this  stage  of 
analysis. 

The  monitored  data  at  sites  12  and  W2  exhibit  no  sig¬ 
nificant  correlation  with  operations  in  their  respective 
vicinities.  In  addition,  they  share  the  same  common 
characteristics  of  weak  or  no  autocorrelation  in  the 
noise  series  and  high  levels  of  variability  relative  lo 
their  mean  levels.  This  raises  a  question  aliout  Ihe 
actual  noise  these  sites  are  measuring;  i.e.,  is  lids 
monitored  noise  rcully  at  rankly  related  to  airport 


KM5.  *1.  Theoretlrnl  nut  no«f  relation  functions. 


operations?  In  terms  of  sampling  requirements,  how¬ 
ever,  these  sites  are  quite  consistent  with  the  other 
sites,  while  the  low  autocorrelation  tends  to  reduce  the 
sample  size  requirements,  the  high  variability  brings 
them  back  in  line  with  those  sites  with  strong  autocor¬ 
relation. 

III.  AN  ALTERNATE  SAMPLING  STRATEGY 

I'ho  rather  large  sampling  requirement  s  revealed  by 
the  analysis  above  suggest  that  an  alternate  strategy 
might  lie  sought  to  reduce  the  amount  of  sampling,. 

When  several  site  s  are  to  be  m<  sutured  it  a  given  fa¬ 
cility,  it  may  become  practical  to  induce  independence 
in  tin-  sample  data  by  spacing  successive  observations 
by  a  sufficient  lag  distance.  The  theoretical  autocor¬ 
relation  function  for  the  ARMA  model  at  a  given  site 
can  be  used  to  estimate  the  spacing  required  to  ap¬ 
proximately  validate  the  independence  assumption. 

For  site  Kl,  the  model  was  AR(1)  with  .ft,  0.5. 

The  theoretical  autocorrelation  lunction  for  an  ARID 

process  is  given  by  p(i  I,  pr  rft*.  A  1.2,3 . 

Hence,  p,  0.5.  p,  0.25,/),  0.  125,  p.,  0.0625.  p, 

0.0312,  etc.,  and  therefore,  tor  practical  purposes, 
observations  spaced  by  5  or  more  lags  are  un  cor  re¬ 
lated.  In  the  case  of  site  A3,  the  model  was  ARMA 
(2, 1)  and  the  theoretical  autocorrelation  function  is 
given  by: 

Po  1  . 

P,  (rt,  -  y.Ml  -  'V  , 

P K  *Np»-»  *  •»  2  . 

Based  on  the  i nixie l  parameters  for  site  Al,  p,  0.28, 
p„  0.23,  p,  0.00,  /»,  0.08,  p,  -0.03  and  again  ob¬ 
servations  spaced  by  5  or  more  days  are  approximated 
uncorrelated.  Figure  4  shows  the  theoretical  autocor¬ 
relation  functions  for  sites  El,  11,  and  W50.  in  de¬ 
veloping  sampling  requirements  based  on  Induced  Inde¬ 
pendence,  the  theoretical  autocorrelations  are  assumed 
to  be  effectively  zero  when  the  correlation  coefficients 
damp  to  less  than  0.05  in  absolute  value. 

A  similar  strategy  can  be  developed  for  randomly 
spaced  groups  of  days  such  as  weekly  blocks  of  data. 

IV.  CONCLUSIONS 

In  the  vicinity  of  airports,  the  data  indicate  that  30 
continuous  days  of  monitoring  is  a  reasonable  estimate 
of  the  number  of  days  required  to  achieve  a  precision 
of  4  2  to  -3  dB  of  the  true  yearly  CNEL  (or  DNL)  value 
with  a  95(?  confidence  level.  Moreover,  the  autocor¬ 
relation  factors  appear  to  be  on  the  order  of  4  to  5  for 
sites  in  the  immediate  vicinity  of  the  airport.  How¬ 
ever,  in  worst -case  situations  such  as  when  the  total 
operations  on  a  runway  become  highly  variable  tsuch 
as  runway  reversals)  when  seasonal  load  or  weather 
factors  become  significant,  or  when  there  is  a  weekly 
cycle  to  the  data,  then  these  numbers  can  become  sig¬ 
nificantly  greater  than  30  and  4,  respectively.  These 
data  indicate  a  worst -rase  requirement  of  60  continuous 
days  in  the  vicinity  of  airports  and  a  worst -case  auto¬ 
correlation  factor  on  the  order  of  8. 
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Because  of  the  autocorrelation  factor  generally  ex¬ 
hibited  in  must  of  the  data  series,  the  number  of  sam¬ 
pling  days  can  be  significantly  reduced  by  inducing  ran¬ 
domness  in  the  selection  of  days  sampled.  That  is, 
sample  days  can  be  selected  sufficiently  far  apart  to 
induce  randomness  in  the  data  gathered,  rather  than 
performing  continuous  monitoring  over  the  total  number 
of  days.  Also,  because  of  the  common  sense  potential 
for  long-term  seasonal  effects,  it  is  recommended  that 
samples  be  selected  from  throughout  the  entire  year. 

A  variety  of  strategies  can  be  employed  based  on  this 
analysis.  For  example,  one  could: 

(a)  Sample  for  a  continuous  period  of  30  days  or  more. 

(b)  Sample  14  days  chosen  randomly  throughout  the 
year,  subject  to  the  constraint  that  no  two  sample  days 
be  less  than  8  calendar  days  apart. 

(c)  Sample  approximately  3  one -week  periods  chosen 
randomly  throughout  the  year,  subject  to  the  constraint 
that  no  two  sample  weeks  be  consecutive. 

The  above  can  be  used  to  achieve  a  precision  of  ♦ 2 
to  -3  dB  of  the  true  yearly  CNEI.  or  DNL  .it  a  95', 
level  of  confidence. 

While  the  above  analysis  pertains  only  to  airports  in 
Southern  California,  subsequent  analysis  currently 
underway  indicate  that  these  sampling  requirements  are 
approximately  valid  In  the  vicinity  of  major  commercial 
airports  on  the  east  coast  as  well. 
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This  paper  is  the  third  in  a  series  dealing  with  the  development  of  temporal  sampling  strategics  for 
estimation  of  mean  noise  levels  in  the  vicinity  of  airports  It  extends  the  previous  analysis  for 
weslcoast,  one-direction  airports  (due  to  prevailing  wmdsl  to  easlcoast,  multidircction  airports 
|  Boston  Logan,  Washington  Dulles,  and  Nalionall  The  results  show  that  the  data  for  many  of  the 
easlcoast  airport  sites  are  nonstationary  in  Ihe  mean  level  and  the  corresponding  tor  iso  uuve 
sampling  requirements  predicted  hy  the  Dynamic  Data  System  |DDS|  methodology  are  veiy 
large,  at  limes  escecdmg  1/1  of  a  year  When  ihe  data  are  stationary.  Monte  Carlo  simulations 
using  the  data  produce  sampling  requirements  comparable  to  the  values  obtained  by  the  DOS 
methodology.  However,  the  DPS  methodology  tends  to  overestimate  sampling  requirements  for 
nonstationary  data  The  simulations  demonstrate  that  nonconseculive  sampling  strategies  reduce 
the  overall  sampling  requirements  for  nonstationary  data.  In  general,  the  results  reveal  the 
following  |a|  Weslcoast  (one-direclionl;  t  50%  precision—  four  weeks,  any  sampling  strategy, 
±  15 %  precision — eight  weeks,  any  sampling  strategy.  |b|  Easlcoast  Imulhdireclionl.  f  b0% 
precision — four  weeks,  one  from  each  quarier,  ^  40%  precision-  -eight  weeks,  one  from  each 
eighth 

PACS  numbers:  43  50.Lj.  43  85. Fm.  4!  50.Qp 


INTRODUCTION 

This  paper  represents  the  third  in  a  series  dealing  with 
temporal  sampling  requirements  for  estimation  for  long¬ 
term  average  sound  levels. The  general  problem  is  the  as¬ 
sociated  statistical  assessment  of  ihe  precision  of  estimates  of 
mean  sound  level.  With  only  a  few  exceptions  such  as  (he 
California  Airport  Noise  Regulation,'  most  techniques  in 
use  today4 '  for  sampling  community  noise  eg  1 1  for  sampling 
over  relatively  short  periods  of  lime,  e  g  .  from  a  few  minutes 
to  perhaps  a  single  day  However,  the  lime  varying  nature  of 
noise  dalu  when  viewed  as  a  lime  series  (hourly  or  dally  aver 
ages/  suggests  I  hat  short  -term  sampling  may  lead  lo  serious 
inaccuracies  m  Ihe  estimation  til  a  long-term  I  yearly)  mean 
noise  level.  Lor  example,  ihe  24-h  periodic  pattern  in  hourly 
mean  sound  level  may  vary  from  about  40  lo  85  dB  The 
Community  Noise  Equivalent  Level  (CNEL)or  Day/Night 
Average  Sound  Level  |DNL)  both  commonly  vary  from  45 
to  80  dB  These  wide  ranges  for  sound  level,  together  with 
the  fact  that  the  data  in  general,  exhibit  high  positive  auto¬ 
correlation  and  high  coefficients  of  variation  suggest  that 
small  and/or  short  sampling  periods  may  provide  both  im 
precise  and  inaccurate  mean  value  estimates 

The  techniques  of  time  series  modeling,  in  general,  pro¬ 
vide  a  powerful  methodology  for  assessment  of  mean  level 
estimation  precision  and  the  formulation  of  sampling  strate¬ 
gies  In  the  first  paper,  the  authors  have  described  the  use  of 
the  Dynamic  Data  System  IDDSI  for  performing  these  anal¬ 
yses  based  on  the  fitting  of  Autoregressive-Moving  Average 
IARMAI  time  series  models  to  the  daily  average  noise  level 
data.  This  previou1  paper  utilized  approximately  one  year's 
daily  CNEL  data  gathered  at  several  sites  around  San  Die- 
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go's  Lindbergh  Field  and  Miramar  Nasal  Air  Station  Anal¬ 
ysis  of  these  data  showed  a  high  degree  of  positive  aulocoi  re¬ 
lation  in  CNEL  values  from  day-to-day.  This  day-to-day 
positive  autocorrelation  is  to  be  expected  because  of  prevail¬ 
ing  winds,  slowly  varying  weather  fronts,  and  the  relatively 
constant  set  of  daily  operations  and  fleet  mix  at  commercial 
airports  This  aulocorrelated  nature  of  the  data,  particularly 


FIG  l  Boston  I  <>gan  International  Airpori  ItX'Jlions  of  noor  monitor 
mg  \»lc* 
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the  degree  of  positive  autocorrelation  between  neighboring 
observations,  increases  the  amount  of  consecutive  sampling 
required  to  estimate  the  long-term  mean  level  with  a  given 
level  of  precision  over  sampling  where  independence  is  as¬ 
sumed 

In  the  second  paper,  the  dynamic  data  system  method 
was  used  to  model  approximately  eight  months  of  daily 
CNEL  data  gathered  at  12  sites  in  the  vicinity  of  the  1-os 
Angeles  International  Airport.  The  results  from  these  first 
two  papers  were  used  to  form  a  set  of  guidelines  for  sampling 
strategies  in  the  vicinity  of  airports.  These  first  papers  used 
C'NEI.  data  and  were  for  weslcoast  airports  only 

Hie  present  paper  extends  the  analysis  to  easlcoast  air¬ 
ports  and  to  the  use  of  ONI..  Specifically,  approximately  1 1 
months  of  daily  ONI.  data  were  obtained  for  1 5  sites  in  the 
vicinity  of  Huston's  Logan  Airport  and  approximately  nine 
months  of  daily  DNLdata  were  obtained  for  nine  sites  in  the 
vicinity  of  Washington's  Dulles  Airport  and  14  sites  in  the 
vicinity  of  Washington's  National  Airport.  The  dynamic- 
data  system  method  was  again  used  to  model  these  daily 
DN  L  data.  Monte  Carlo  simulations  were  performed  to  ver¬ 
ify  the  sampling  requirements  obtained  from  the  DOS  mod¬ 
eling  of  the  data  and  to  study  alternatives  to  consecutive 
sampling.  The  results  of  these  analyses,  along  with  the  re¬ 
sults  of  the  two  previous  papers  are  used  to  form  a  set  of 
guidelines  for  sampling  strategies  in  the  vicinity  of  airports. 


FIG  2  Washington  National  Airport  —  knatiom  of  mwv  monitoring  otn 
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FIG  .3.  Washington  Dulles  Airpori  —locations  of  noise  momionng  sites 


I.  THE  DATA  BASE 

Continuous  daily  monitored  DNL.  values  were  supplied 
by  Logan  Airport  for  the  period  from  I  October  1478  to  31 
October  1979  for  the  15  monitoring  locations  indicated  in 
Fig.  I.  These  daily  DNL  data  were  transformed  to  values 
denoted  by  X,  on  a  linear  scale,  proportional  to  time-weight¬ 
ed  daily  sound  exposure  as  shown  in  Eq  1 1 1, 

x  —  i(y,>NL/l0'  ,i, 

As  noted  in  earlier  papers,  these  X  values  arc  used  because 
one  is  interested  in  estimating  a  yearly  mean  DNL  or  CNEL 
value  which  by  definition  is  estimated  by  the  sample  yearly 
mean  (arithmetic  average)  X  Monitored  data  are  frequently 
measured  for  far  less  than  one  year's  time  period,  and  it  is  the 
purpose  of  this  research  to  address  the  validity  of  techniques 
to  estimate  this  yearly  mean 

The  time  senes  of  X  values  at  some  of  the  locations  were 
divided  into  part  A  and  part  B  at  a  gap  in  the  data  around  the 
halfway  mark.  Due  to  sizeable  gaps  in  some  of  the  sites'  data, 
only  about  one-half  year  of  continuous  data  were  available 
and  modeled.  Using  the  DDS  method,  ARMA  models  of 
vanous  orders  were  developed  for  15  of  the  16  above  time 
senes  of  the  quantity  X  (For  a  bnef  review  of  the  modeling 
met  hods  and  terminology,  the  reader  is  referred  to  Appendix 
A  )  The  estimated  parameters  were  used  to  determine  the 
autocorrelation  factor,  coefficient  of  variation,  and  appro¬ 
priate  sampling  requirements  for  mean-level  estimation 
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The  partitioning  of  the  data,  as  mentioned  above,  pro¬ 
vides  for  the  opportunity  to  examine  the  homogeneity  of  the 
stochastic  structure  of  the  data  over  the  entire  period  of  one 
year  Although  the  DDS  models  do  not  reveal  seasonal  var¬ 
iations.  modeling  data  in  different  times  of  the  year  does 
show  some  differences  in  stochastic  structure  at  certain  sites 
This  fact  leads  to  the  postulation  of  several  alternate  sam¬ 
pling  strategies  including  sampling  in  each  of  the  four  sea¬ 
sons.  This  will  be  discussed  in  detail  later  in  the  paper 

Runway  operations  data  were  unavailable  at  Boston 
Logan,  but  all  the  sites  were  evaluated  and  the  monitored 
levels  found  to  exceed  reasonable  estimates  of  the  communi¬ 
ty  noise  in  the  absence  of  aircraft  Thus  the  monitored  data 
were  assumed  to  be  predominantly  aircraft  noise  and  subject 
to  a  constant  set  of  prevailing  operations 

Data  were  analyzed  from  Washington  National  and 
Dulles  Airports  for  the  penud  of  17  March  W7H  to  14  De¬ 
cember  W78.  At  Dulles  Airport,  the  daily  number  of  arrivals 
and  departures  are  significantly  less  than  at  Los  Angeles  nr 
Boston,  and  for  both  Washington  Airports,  the  monitoring 
sites  were  generally  positioned  farther  from  the  airport  (see 
Figs  2  and  3).  For  the  Washington  Airports,  it  was  attempt¬ 
ed  to  determine  whether  each  site  was  monitoring  predomin¬ 
ately  aircraft  or  community  noise.  Monthly  Federal  Avi¬ 
ation  Administration  noise  level  plots  (Fig.  4)  served  as  the 
basis  for  this  classification  process.  The  sample  plot  in  Fig.  4 
shows  the  airport  noise  level  (indicated  by  the  LEQA  hatch¬ 
ing)  standing  clearly  above  the  background  community 
noise  level  of  about  6$  dB.  Plots  of  this  type  were  studted  for 
all  of  the  National  and  Dulles  sites.  In  some  cases,  a  site 
could  be  unquestionably  classified  as  either  an  airport  or  a 
community  noise  site,  but  in  other  cases  it  was  necessary  to 
classify  a  site  as  a  mixed  airport  and  community  noise  site. 

II.  DISCUSSION  OF  MODELING  RESULTS 
A.  Boston's  Logan  aftss 

For  Boston's  Logan  Airport,  frequency  histograms  and 
lime  senes  plots  were  developed  and  can  be  found  in  Ref.  6. 


The  DDS  modeling  results  are  listed  in  l  able  I  for  all  sites  .u 
Boston  For  those  sites  whose  data  were  divided  inloappin 
imately  two  equal  parts,  modeling  results  are  given  for  ea.  :t 
part  and  for  the  entire  data  sel  Mosl  of  the  dala  were  I  mm  ,.1 
to  follow  ARllltime  series  models,  w  ith  autocorrelation  t.u 
tors  commonly  ranging  from  I  t  to  3  0  The  required  samp.e 
sizes  for  estimation  of  the  mean  within  a  ,  fO't  ol 
Jtl  +  2,  -  3  dB|  at  the  95r'<  confidence  level  commonly  var 
led  from  10  to  XO  days  Sne4A  was  the  most  extreme  excep 
lion,  with  a  very  high  autocorrelation  factor  of  6  '5  and  a 
correspondingly  large  sample  size  of  1 14  days  Sites  1  and  4 
fin  the  same  general  areal  also  have  higher  than  average  auto¬ 
correlation  factors  Almost  all  Ihe  coefficients  of  variation 
range  from  0  6  to  I  3 

The  dala  which  were  divided  into  two  pans  arc  used  m 
more  carefully  examine  the  long-term  i  yearly  i  sound  charac 
lerisltcsl  I  able  1 1  In  some  cases  Isites  1 .  V  and  Hi,  the  model 
mg  results  arc  very  similar  in  terms  of  both  autocorrelation 
factor  and  coefficient  of  variation,  going  rise  to  very  similar 
sampling  requirements  For  others  isites  t.  4.  and  hi.  the 
results  are  vimewhat  different  Of  particular  interest  is  ihe 
fact  that  modeling  of  the  combined  data  for  sites  3.  4.  and  6 
produces  sampling  requirements  generally  larger  than  for 
either  part  For  site  3.  the  main  difference  in  the  modeling 
results  is  the  autocorrelation  factor  For  this  site,  exatmna 
lion  of  the  data  (Fig.  5(a|]  reveals  a  gradual  upward  trend 
from  the  middle  of  the  second  part  of  the  dala  This  nonsta¬ 
tionary  tendency  increases  Ihe  order  of  the  model  and  pro¬ 
duced  a  higher  autocorrelation  factor  due  to  a  model  root 
approaching  1 .0.  For  site  6.  the  second  half  of  the  data  has  a 
more  random  structure  with  corresponding  autocorrelation 
factor  equal  to  10  In  examining  site  4.  it  is  seen  in  Table  I 
that  the  model  for  the  first  part  of  the  year  gives  nse  to  a  very 
large  sampling  requirement,  134  days.  This  is  2.5  times 
greater  than  any  other  site  Examination  of  the  data  [Fig 
5(b|)  shows  a  large  downward  shift  in  the  mean  level  of  thc 
data  after  about  the  first  two  months  This  nonstationarity 
has  produced  a  model  with  high  autocorrelation  factor  Al- 


FIG  4  Typical  FA  A  monthly  data 
summary  at  a  Washington  monitonng 
sue  In  this  FAA  figure,  ’  LEQA"  are 
the  equivalent  levels  (by  hour!  resulting 
from  aircraft  noises  ”LEQ<  "  are  the 
equivalent  levels  resulting  from  I  nonair 
craft)  community  noises,  and  LEO  are 
the  total  equivalent  levels  resulting  from 
the  combination  of  community  and  air¬ 
craft  noises 
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though  no  records  art-  available  to  (he  authors  such  a  data 
(rend  could  he  caused  by  a  marked  and  sustained  change  in 
runway  activity  or  flight  patterns  or.  perhaps  instrumenta¬ 
tion  problems.  The  combined  model  appears  dominated  by 
the  trend  in  the  first  part  of  the  data  and  shows  large  sam¬ 
pling  requirements. 

B.  Washington  Dulles  and  National  sitaa 

DDS  modeling  was  performed  for  each  site  in  the  vicin¬ 
ity  of  Dulles  and  National  Airports.  Site  classification  and 
modeling  results  are  summarized  in  Tables  II.  Ill,  and  IV 
lairport  sites,  community  sites,  and  mixed  sites,  respective¬ 
ly)  Low-order  models  resulted  for  most  of  the  sites,  with  a 
first-order  autoregressive  model,  AR(  1 1,  being  common.  The 
sampling  requirements  for  all  the  Washington  sites  are  gen¬ 
erally  similar  to  those  for  the  Logan  sites. 

Tables  II,  III,  and  IV  indicate  that  several  Washington 
Dulles  and  National  sites  exhibit  wide  differences  in  model¬ 
ing  results  when  each  part  and  the  combined  data  are  exam¬ 
ined  separately.  In  particular,  for  seven  sites  (Dulles  I,  6,  7, 
and  lOand  National  14,  21,  and  22).  DDS  modeling  results 
vary  considerably  from  one  part  of  the  data  to  the  other  and 
sampling  requirements  based  on  the  combined  data  models 
arc  generally  larger  than  for  either  part  Examination  of  the 
lime  series  data  show  nonstationary  (rends  either  within  or 
across  both  parts  of  the  data.  These  trends  can  again  he  char¬ 
acterized  as  a  shifting  of  the  mean  level  of  (he  data.  The 
extreme  case  is  Dulles  sue  h  which  lias  small  sampling  re¬ 
quirements  for  each  half  (A’  -  5|  bill  a  very  large  require¬ 
ment  l,V  1 58|  for  the  combined  data  I  nix  sue  s  data  (Fig 
S|c||  show  a  marked  shift  in  the  mean  level  at  about  (lie  mid¬ 
year  point 
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III.  COMPARISON  OF  EASTCOAST  AND  WESTCOAST 
AIRPORTS 

In  summarizing  the  modeling  resulls  across  ihc  liostou 
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I  AM  t  II  Washington  National  and  Dulles  Airport  sues  — DOS  modeling  results  on  original  data 


Airport 

Site 

Model 

Coefficient 

of 

variation 

Autocor  relation 
factor 

Sample 

sue 

Dulles 

7A 

ARID 

051 

1  49 

7 

Dulles 

7B 

AR|I| 

0  65 

2  50 

P 

Dulles 

7 

ARMA|3.  )) 

0.58 

6  17 

34 

National 

1  3A 

White  noiae 

1  83 

1  00 

54 

National 

13B 

AR|2) 

066 

3  25 

23 

National 

13 

ARMAI4.  3| 

1  40 

1.47 

48 

National 

I4A 

AR|1| 

047 

1  49 

6 

National 

MB 

ARID 

0  62 

309 

19 

National 

M 

ARMAI2.  2) 

061 

16  38 

97 

National 

I5A 

While  noise 

0.86 

1  00 

13 

National 

I3B 

ARID 

0  80 

l  79 

19 

National 

15 

arid 

083 

1  36 

15 

National 

16 

ARID 

0  53 

2  56 

12 

National 

19 

ARID 

0  92 

1  82 

25 

fABLE  Ill  Washington  NatMHiai  and  Dulles  community  sites — DDS  modeling  results. 


Airport 

Site 

Model 

Coeffictenl 

of 

vanation 

Autocorrelation 

factor 

Sample 

size 

Dulles 

IA 

While  none 

032 

too 

5 

Dulles 

IB 

ARID 

0.95 

206 

31 

Dulles 

l 

ARMAll.  I| 

088 

7.24 

90 

Dulles 

4 

Nonstationary 

0  53 

Dulles 

5 

AR(3| 

064 

1028 

67 

Dulles 

9 

ARID 

048 

1.77 

7 

National 

24 

ARID 

043 

5.14 

14 

TABLE  IV  Washington  National  and  Dulles  mixed  sites — DDS  modeling  results 


Airport 

Site 

Model 

Coefficient 

of 

vanation 

Autocorrelation 

factor 

Sample 

size 

National 

3B 

AR|2) 

0.36 

2.65 

6 

Dulles 

6A 

White  noise 

1  47 

100 

5 

Dulles 

6B 

ARID 

037 

2.11 

5 

Dulles 

6 

ARMAll.  1) 

1  44 

4  78 

158 

Dulles 

8A 

White  noise 

1.55 

1  00 

9 

Dulles 

8B 

White  notse 

2  02 

1  00 

II 

Dulles 

8 

White  noise 

071 

1  18 

10 

Dulles 

10A 

ARID 

0  76 

2  62 

25 

Dulles 

106 

Nonstaimnary 

062 

Dulles 

10 

ARMAI3.  3) 

0  70 

8  02 

62 

National 

II 

ARMAI2,  2) 

080 

n  no 

112 

National 

12 

AR|2| 

0  65 

3  80 

26 

National 

18 

White  noise 

0  92 

1  00 

14 

Nai  tonal 

20A 

While  noise 

0  83 

1  00 

11 

National 

206 

ARID 

051 

1  67 

7 

National 

20 

White  noise 

0  74 

1  os 

10 

National 

2IA 

Nonstationary 

0  62 

National 

2IB 

ARID 

103 

1  85 

32 

National 

21 

ARMAI3,  1) 

0.62 

5  51 

67 

National 

22A 

ARID 

0  83 

1  30 

15 

National 

22B 

ARI2I 

080 

483 

50 

National 

22 

ARMAll.  1) 

082 

14  79 

159 

s 

I 
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beroflhe  sites  exhibit  nonstationary  behavior,  viz.,  changing 
mean  level  over  the  year.  As  a  result,  long-term  consecutive 
sampling  requirements  are  very  large,  often  constituting 
more  than  one-third  of  a  year.  This  result  is  much  different 
than  the  sampling  requirements  analysis  for  the  westcoast 
airports  which  exhibit,  in  general,  stationary  stochastic 
structure  over  an  entire  year’s  data.  In  attempting  to  delin¬ 
eate  differences  in  the  characteristics  of  eastcoast  and  west- 
coast  airports  and  the  sampled  data  obtained,  two  observa¬ 
tions  can  be  noted:  (1)  the  westcoast  airports  are 
one-runway-direction  airports,  and  (2)  the  monitoring  sites 
for  the  westcoast  airports  are  generally  located  closer  to  the 
runways. 

The  results  relative  to  the  eastcoast  airports  suggest 
that  any  analysis  should  be  confined  to  data  covering  sub¬ 
stantially  a  period  of  a  year.  Thus  the  remainder  of  this  paper 
will  focus  only  on  sites  for  which  a  full  year's  data  were 
available. 

Figure  6  summarizes  the  results  for  all  of  the  airports 
modeled  (airport  sites  only),  including  Dulles,  National,  and 
Logan  from  this  paper  and  Los  Angeles  (including  one  site 
from  Lindberg  Field)  from  the  author's  previous  papers. 
This  figure  was  constructed  to  more  graphically  represent 
the  similarities  and  differences  among  the  airports  in  terms 
of  their  sampling  characteristics.  The  results  generally  show 
that  the  westcoast  airports  (typically  one-direction,  due  to 
prevailing  winds  off  the  ocean)  lend  to  have  lower  coeffi¬ 
cients  of  variation  and  comparable  autocorrelation  factors 
relative  to  the  multirunway  and/or  multidirection  (variable 
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AUTOCORRELATION  FACTORS 
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FIO  6.  Airport  modeling  results.  The  “D"  in  the  Washington  Airport  data 
indicates  Dulles  tiles  and  the  '*S"  in  the  Los  Angeles  Airport  data  indicates 
the  Lindberg  Field  site. 


witnl)  raslconst  jtrporls  These  results  produce  ovettilt  sum 
plmg  requirements  for  the  weslcoast  airports  which  are  gen 
erally  lower  than  those  for  the  eastcoast  utrports 

Data  and  modeling  results  from  the  Dulles  anti  Nation 
al  community  sites  and  mixed  sites  arc  surprising  in  that  they 
do  not  appear  to  differ  significantly  from  ihe  airport  data 
Thai  is.  the  autocorrelation  factors,  coefficients  of  variation 
and  therefore  overall  sample  size  requirements  do  not  ap¬ 
pear  to  differ  significantly  among  types  of  sites 

IV.  MONTE  CARLO  STUDY  OF  GENERALIZED 
SAMPLING  STRATEGIES 

Because  of  the  presence  of  nonslat tonary  trends  and 
large  sampling  requirements  for  some  of  the  data,  Monte 
Carlo  sampling  experiments  were  performed  with  Ihe  Los 
Angeles,  Boston,  and  Washington  data  Through  such  simu¬ 
lations.  generalized  sampling  strategies  including  alterna¬ 
tives  to  consecutive  sampling  may  be  examined  Such  alter¬ 
nate  strategies  may  require  fewer  total  samples  than 
consecutive  sampling  and  provide  a  means  to  accommodate 
trends  in  the  airport  noise  data  over  a  period  of  about  one 
year. 

Sampling  experiments  were  performed  on  those  sues 
with  one  year  of  reasonably  continuous  data  In  ihe  first  sel 
of  experiments,  Ihe  total  number  of  samples  taken  to  csti 
male  the  mean  noise  level  for  each  strategy  was  2X  days 
Four  sampling  strategies  were  investigated:  (1)28  days  ran¬ 
domly  spaced  throughout  the  year.  (2)  one  week  of  consecu¬ 
tive  sampling  in  each  quarter  of  the  year,  (3|  two  weeks  of 
consecutive  sampling  in  each  half  of  the  year,  and  |4|  28 
consecutive  days. 

For  each  strategy,  the  starting  date  was  selected  ran¬ 
domly  for  each  sample  or  group  or samples  For  each  site,  the 
sampling  strategy  was  repeated  20  times  For  each  Inal,  the 
sample  mean  noise  level  was  calculated,  and  the  variability 
in  means  among  the  trials  was  used  to  estimate  the  variance 
of  sample  means.  This  variance  estimate  together  with  the 
appropriate  f  statistic  (/,„„,„  =  2.039)  was  used  to  develop 
an  estimate  of  predictive  precision  as  a  percentage  of  the 
population  mean  at  the  a  =  0.03  significance  level.  For  the 
time  series  approach,  it  was  desired  to  determine  the  require¬ 
ments  for  consecutive  sampling  to  estimate  the  mean  noise 
level  within  ±  50%  of  the  mean  at  the  95%  confidence  lev¬ 
el.  For  these  Monte  Carlo  sampling  experiments,  the  percent 
precision  P  for  each  sampling  strategy  is  computed 

In  a  second  set  of  experiments  a  total  of  56  samples  were 
taken.  Five  sampling  strategies  were  considered:  ( 1 1  56  days 
randomly  spaced  throughout  the  year,  (2)  one  week  of  con¬ 
secutive  sampling  in  each  eighth  of  the  year,  (3)  two  weeks  of 
consecutive  sampling  in  each  quarter  of  the  year,  (4)  four 
weeks  of  consecutive  sampling  in  each  half  of  the  year,  and 
(5)  56  consecutive  days. 

Figure  7(a)  through  7(g)  show  the  Monte  Carlo  simula¬ 
tion  results  for  Los  Angeles,  Boston,  National,  and  Dulles 
Airports,  respectively.  In  examining  these  figures,  a  number 
of  important  observations  are  noted.  For  Los  Angeles,  peri¬ 
odic  sampling  does  indicate  a  slight  but  not  marked  im¬ 
provement  in  predictive  precision  over  consecutive  sam¬ 
pling.  With  the  exception  of  two  sites,  the  Los  Angeles 
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result*  show  that  a  ±  50%  precision  can  be  attained  with  28 
samples,  regardless  of  the  sampling  strategy  chosen.  For 
Boston  Logan  and  the  Washington  Airports  signiAcanl  im¬ 
provements  in  the  predictive  precision  can  be  achieved  by 
periodic  sampling,  e.g.,  one  week  from  each  quarter  over  the 
year.  This  is  particularly  true  for  those  sites  which  exhibited 
nonstationary  behavior.  To  guarantee  a  ±  60%  precision 
level  for  these  airports  it  is  required  to  sample  one  week  from 
each  quarter  over  the  entire  year. 

In  considering  the  results  of  the  simulation  experiments 
involving  requirements  of  36  days  of  sampling  it  is  noted  that 
for  Los  Angeles,  ±  35%  predictive  precision  is  obtained  for 
all  sites  regardless  of  the  sampling  strategy  except  for  sites  12 
and  W3.  For  these  two  sites  ±  33%  precision  is  attainable 
by  sampling  for  one  week  out  of  each  eighth  of  the  year  For 
the  Boston  and  Washington  Airports,  ±  40%  precision  can 
be  achieved  for  all  sites  if  eight  one-week  samples  are  taken, 
one  from  each  eighth  of  the  year. 

For  Lot  Angeles,  the  DDS  modeling  consecutive  sam¬ 
pling  requirements  and  those  obtained  from  the  Monte 
Carlo  simulations  are  generally  about  the  same  (see  Table  V). 
These  sites  exhibit  stationary  stochastic  structure  for  the  en¬ 
tire  year  and  the  simulations  verify  the  credibility  of  the 
DDS  modeling  results.  The  same  comparison  holds  for  Bos¬ 
ton  sites  1,  5,  8,  Dulles  site  8,  and  National  sites  15  and  20. 
These  sites  exhibit  stationary  behavior. 


TABLE  V  Companion  of  DDS  and  Monte  Carlo  simulation  results 


Site 

DDS  modeling 
results 

consecutive 
samples  for 

P-  ±50% 

Monte  Carlo 
simulation 
results  %P  for  28 
consecutive 
samples 

Monte  Carlo 
simulation 
results  %P  for  56 
consecutive 
samples 

LAX 

At 

<> 

260 

17.0 

LAX 

A  2 

19 

32.0 

220 

LAX 

El 

4 

240 

13  0 

LAX 

E3 

8 

320 

24.0 

LAX 

12 

30 

770 

470 

LAX 

L2 

II 

350 

180 

LAX 

W2 

16 

360 

25  0 

LAX 

W3 

60 

81.0 

320 

LAX 

W4 

II 

400 

280 

Boston 

1 

II 

400 

230 

Bow  on 

3 

60 

710 

780 

Bom  on 

4 

128 

1210 

1 160 

Bom  on 

5 

35 

63.0 

490 

Bom  on 

6 

79 

93.0 

420 

BoMon 

8 

31 

510 

390 

Dulles 

1 

90 

37.0 

460 

Dulles 

4, 

Nonstationary  74.0 

47.0 

Dulles 

6 

158 

93.0 

56.0 

Dulles 

7 

34 

470 

240 

Dulles 

8 

10 

32.0 

20.0 

Dulles 

10 

62 

53.0 

390 

National 

13 

48 

720 

460 

National 

14 

97 

55.0 

320 

National 

IS 

15 

47.0 

350 

National 

18 

14 

960 

900 

National 

20 

10 

260 

20  0 

National 

21 

67 

740 

570 

National 

22 

159 

800 

660 
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For  the  bites  exhibiting  nonslationary  behavior  at  Bos 
ton,  Dulles,  and  National,  the  comparison  of  the  DDS  mod 
cling  and  simulation  results  are  not  always  consistent  In 
particular,  the  DDS  models  may  overquote  the  consecutive 
sampling  requirements  necessary  to  achieve  a  particular  lev  - 
el  of  precision  when  compared  with  the  simulation  results 
This  is  particularly  true  for  site  6  at  Boston,  site  6  at  Dulles 
and  sites  14  and  22  at  National  For  a  comprehensive  exami¬ 
nation  of  the  Monte  Carlo  simulation  results,  the  reader  is 
directed  to  the  Appendix. 

V.  CONCLUSIONS 

The  results  generally  show  that  the  westcoast  airports 
{typically  one-direction  due  to  prevailing  winds)  tend  to  have 
lower  coefficients  of  variation  and  comparable  autocorrela¬ 
tion  factors  relative  to  the  multirunway  and/or  multidirec 
lion  {variable  wind)  eastcoast  airports  These  results  produce 
overall  consecutive  sampling  requirement  lor  the  westcoast 
airports  which  are  generally  lower  than  those  for  the  east 
coast  airports. 

Many  of  the  eastcoast  airport  sites  exhibit  nonstation- 
ary  trends.  This  nonstationary  condition  is  evidenced  by 
vastly  different  models  and  resulting  sampling  requirements 
for  the  data  as  whole  or  for  one  part  of  the  data  when  com¬ 
pared  with  the  other  part.  Typically,  the  sampling  require¬ 
ments  for  the  entire  year's  data  greatly  exceed  the  require¬ 
ments  derived  for  one  or  both  parts 

Monte  Carlo  simulations  using  the  data  show  similar 
results  to  the  DDS  methodology  when  the  data  are  station¬ 
ary.  However,  the  DDS  methodology  overestimates  the 
sampling  requirements  for  nonstationary  data  Moreover, 
the  Monte  Carlo  simulations  show  that  nonconsecutive  sam¬ 
pling  strategies  for  nonstationary  data  reduce  the  overall 
sampling  requirements. 

These  results  can  be  generalized  to  the  eastcoast  and 
westcoast  airports  as  follows:  (a)  Westcoast  lone-direction); 
±  50%  precision — four  weeks,  any  sampling  strategy. 
±  35%  precision — eight  weeks,  any  sampling  strategy.  |b) 
Eastcoast  (multidirection);  ±  60%  precision — four  weeks, 
one  from  each  quarter.  ±  40%  precision — eight  weeks,  one 
from  each  eighth. 

In  general,  the  stationary  airport  sites  are  modeled  as 
AR(I)  processes.  Occasionally,  there  are.  white  noise  sites 
and,  occasionally,  there  are  nonstationary  sites.  The  com¬ 
munity  sites  and  the  mixed  sites  also  are  generally  AR(I) 
models. 
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APPENDIX  A:  STOCHASTIC  MODELING  BY  THE 
DYNAMIC  DATA  SYSTEM  (DDS)  METHOD 

The  method  of  dynamic  data  system  IDDS)  provides  for 
the  development  of  parametric  stochastic  time  senes  models 
of  the  autoregressive — moving  average  IARMAI  class  Dai¬ 
ly  sound  exposure  [Eq  (1)]  when  viewed  as  a  time  senes  of 
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values  *,,  X2,...JCN  has  been  shown  to  be  well  characterized 
by  such  models''2  which  makes  it  possible  to  determine  the 
precision  associated  with  an  estimate  of  the  yearly  mean 
sound  exposure  level  when  the  observed  daily  values  X,  are 
autocorrelated 

The  general  ARM  A  model  is  given  by 

X,  =  4,X,  i  +  4iX,  »+■■•  +  4. X,  .  +  a, 

-  0,o,  ,-0-fi,  1 - 0.0,  „  . 

where  X,  is  the  noise  level  (daily  average)  for  day  l,  a,  is  the 

random  disturbance  for  day  /.  4 . 4.  are  autoregressive 

parameters,  and,  tf . ,0m  are  moving  average  parameters. 

Given  the  time  series  X,,  the  appropriate  order  of  the 
AKMA  model,  viz.,  the  proper  values  for  n  and  m,  may  he 
determined  and  the  parameters  of  the  model  may  be  estimat¬ 
ed  by  ihe  method  of  least-squares  f:or  details  of  the  modeling 
procedures,  the  reader  is  referred  to  Ref.  7.  The  majority  of 
the  filled  ARM  A  models  obtained  for  the  data  analyzed  in 
this  paper  are  of  relatively  low  order,  e  g.,  , 


AR(1).  X,  =4,X,  ,+fl,  , 

AR(2|:  X,  =  4, X,  .  ,  +  4J, 

ARM A(  1,1):  X,  =  4,X,  ,  +  a,  -B,a,  , 

The  ARMA  models  fit  to  the  daily  sound  exposure  time 
series  X,  may  be  used  to  estimate  the  precision  associated 
with  the  sample  mean  X.  It  can  be  shown1  that  the  variance 
estimate  of  the  sample  mean  is  given  by 

Variance**)  -  ^-  [(>  -  £  *,)/(>  -£<*,)]  . 

where  a3,  is  the  residual  mean  square  for  the  fitted  ARMA 
model.  Given  the  above  variance  estimate,  100(1  -  a)% 
confidence  intervals  of  the  form 

J  ±  'v  ,  m,  I  „/!  I  Variance!*  | ) 1  /! 
may  be  obtained  for  the  true  yearly  mean  sound  exposure 
level,  thereby  providing  anestimate  of  the  precision  associat¬ 
ed  with  the  sample  mean  * 
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Site 

One  week 

four  times  per  year 

%P 

Sampling  strategy 

Two  weeks 
two  times  per  year 
%/> 

28 

consecutive 

days 

%P 

28 

random 

days 

%P 

Los  Angeles 

A1 

20.0 

23.0 

260 

15  0 

Los  Angeles 

A2 

21.0 

20.0 

32.0 

14.0 

Los  Angeles 

El 

14.0 

18.0 

24.0 

11.0 

Los  Angeles 

E2 

13.0 

22.0 

32.0 

10.0 

Los  Angeles 

II 

44.0 

54.0 

490 

270 

Los  Angeles 

12 

52.0 

83.0 

77.0 

49.0 

Los  Angeles 

L! 

16.0 

25.0 

37.0 

12.0 

l.iis  Angeles 

L2 

31.0 

35.0 

35.0 

20.0 

Los  Angeles 

W2 

48.0 

42.0 

360 

58.0 

Los  Angeles 

W3 

78.0 

81.0 

81.0 

45.0 

Los  Angeles 

W4 

24.0 

46.0 

40.0 

210 

Boston 

1 

24.0 

37.0 

40.0 

21.0 

Boston 

3 

40.0 

66.0 

71.0 

39.0 

Boston 

4 

63.0 

96.0 

121.0 

59.0 

Boston 

5 

54.0 

59.0 

65.0 

39.0 

Boston 

6 

53.0 

77.0 

93.0 

57.0 

Boston 

8 

56.0 

43.0 

510 

38.0 

Dulles 

1 

41,0 

57.0 

57.0 

300 

Dulles 

4 

44.0 

62.0 

74.0 

16.0 

Dulles 

6 

55.0 

64.0 

93.0 

30.0 

Dulles 

7 

32.0 

37.0 

47.0 

21.0 

Dulles 

8 

35.0 

32.0 

32.0 

28.0 

Dulles 

10 

39.0 

55.0 

550 

25.0 

National 

13 

62.0 

74.0 

72.0 

570 

National 

14 

25.0 

50.0 

55.0 

21  0 

National 

IS 

29  0 

45.0 

47  0 

30  0 

National 

IK 

*4.0 

94.0 

960 

440 

N  .t:>':  :) 

.!() 

:>5  o 

29.0 

26  0 

4.3.0 

Vll'OI  .1 

.'1 

MO 

t>4  0 

74  0 

25  0 

NaltiM.il 

wn 

MM) 

KOI) 

2.3.0 

SM:> 


J  AtouM  i»fH  Mm  Vol 


N»l  (>.  e/wm*  1 963 


Schom#r  §1  Sampling  sir ilexes -  airports 


2049 


Sampling  strategy 


Sue 

One  week 
eight  lime* 
per  year 

%P 

Two  weeks 
four  times 

per  year 

%P 

f  our  weeks 

two  times 

per  year 
%P 

56 

consecutive 
days 
<yt  p 

56 

random 

days 

7rP 

Los  Angeles 

A  1 

12.0 

19  0 

20.0 

17  0 

10  0 

1  os  Angeles 

A2 

1  1  0 

24.0 

21.0 

22  0 

!  1.0 

I  on  Angeles 

LI 

II  0 

1 10 

15.0 

15  0 

:oo 

l.os  Angeles 

L2 

12(1 

160 

21  0 

24.0 

90 

1  ns  Angeles 

II 

24  (1 

.lo.o 

17.0 

10  0 

25  0 

1  ns  Angeles 

1 2 

14  II 

15.0 

2>M) 

47  0 

15.0 

1  os  Angeles 

1  1 

Kill 

12  0 

210 

1ft  0 

8  O 

Lift  Angeles 

1.2 

18.0 

20  0 

19.0 

18  0 

17.0 

l  os  Angeles 

W2 

24.0 

28  0 

29,0 

25.0 

19  0 

.  .os  Angeles 

W1 

11  0 

190 

47  0 

52  0 

lft.O 

i.ns  Angeles 

W4 

18.0 

290 

300 

28.0 

1  7.0 

ikiston 

1 

21.0 

15  0 

25.0 

250 

160 

Boston 

.1 

500 

15.0 

18.0 

78.0 

16.0 

Boston 

4 

19.0 

47.0 

56.0 

1 160 

25  0 

Boston 

5 

.14.0 

27.0 

27.0 

49.0 

25  0 

Boston 

0 

27.0 

41.0 

19.0 

420 

11.0 

Boston 

8 

260 

15.0 

25.0 

19  0 

21.0 

Dulles 

1 

27  0 

.10.0 

40.0 

46.0 

29.0 

Dulles 

4 

2.1 0 

14.0 

47.0 

47.0 

14.0 

Dulles 

0 

.180 

51.0 

290 

560 

14.0 

Dulles 

7 

160 

21.0 

21.0 

240 

15.0 

Dulles 

X 

17.0 

11.0 

27.0 

20.0 

22.0 

Dulles 

III 

270 

.12.0 

50.0 

29.0 

lft.0 

National 

1.1 

12.0 

60.0 

410 

46.0 

40.0 

National 

14 

17  0 

25.0 

19.0 

520 

19.0 

National 

15 

19.0 

23.0 

12.0 

15.0 

24.0 

National 

IK 

15  0 

56.0 

750 

90.0 

32.0 

National 

20 

20.0 

21.0 

21.0 

200 

22.0 

National 

21 

28.0 

58.0 

58.0 

570 

24.0 

National 

22 

18.0 

40.0 

56.0 

66.0 

18.0 

’R  E.  DeVor,  P  D  Schomet,  W  A  Kline,  ind  R  D.  Ncalhamer.  “Devel¬ 
opment  of  Temporal  Sampling  Strategies  for  Monitoring  Noise,"  J 
Acoust  Soc  Am.  M.  763-771  (1979). 

:P  D  Schomer  and  R.  E.  DeVor.  “Temporal  Sampling  Requirements  for 
EatimatKin  of  Long-Term  Average  Sound  Levels  in  ihc  Vicinity  of  Air- 
porta,"  J.  Acoust  Soc.  Am.  713-719(1981). 

'State  of  California,  Department  of  Aeronautic*.  Title  21.  Subchapter  b, 
None  Standards,  28  November  1970.  Rev  22  June  1979 

‘D  C.  Pie*  and  L.  C  Sutherland.  "Evaluation  of  Spatial  Sampling  Tech¬ 
nique*  for  Community  Notse  Surveys,"  Wyle  Research.  El  Segundo.  CA, 


Rep  WR  77-5.  U  S.  EPA  (1977) 

'J.  Stearns  et  a! ,  "Community  Noise  Monitoring— A  Manual  for  Imple¬ 
mentation,"  Wyle  Research,  El  Segundo.  CA,  Rep  WR  76- 8,  U  S  EPA 
(19761 

*P  D  Schomer  et  al ,  "Temporal  Sampling  Requirement*  for  Estimating 
the  Mean  Level  in  the  Vicinity  of  Military  Installations."  U  S  Army 
Const  Eng.  Res.  Lab.  Tech.  Rep.  N-101  (April  1981). 

'S  M  Pandit  andS.  M.  Wu,  Time  Series  and  System  Analysts  Modeling  and 
Applications  (Wiley ,  New  York,  1982) 


2050 


J  Acoust  Soc  Am  .  Vd  73.  No  6.  Juno  1 963 


Schomer  or  a /  Sampling  strategies— airports 


2050 


I  Ml.  UlsrMtttlil  I'KM 


I  Met  uf  I  nyllMMI  I 
AriNi  I*  It  MlMt  tu« 

MIMi  IMf M-ASl-l.  I?) 

AXIHt  DA£H-CCP 
Arm  i  daen-cw 
Arm  i  dacn-cwe 
ATm  i  DACN-CWH-R 
Arm i  otfM-cw) 

ATm  I  OACN-CWP 
ATTN.  DACN-EC 
ATm  I  DAEH-ECC 
ATm  I  DACN-ECE 
ATm»  cw«H-zrf 
ATm  i  0A£N-€C8 
ATm  I  OAEN-RD 
ATm  I  OACN-RDC 

Arm  i  D^N-RDN 
Arm  I  dacm-rn 
ATm.  0ACN-2C2 
Arm*  o**-zce 
ATm  I  OA£N-ZCI 
ATm*  OAEN-ZCM 

FESA,  ATm  I  Library  22060 
Arm.  OCT  III  79906 

US  A my  Engineer  Districts 
ATm.  Library  <  41  > 

US  Amy  Engineer  Divisions 
ATm.  Library  (14) 

US  Amy  Europo 

AfcA£N-QOC$/Engr  09403 
I SAC  09061 
V  Corps 

Aim.  och  (in 
VII  Corps 
ATm.  OCH  (IS) 

21st  Support  Comma d 
ATm.  OCH  (12) 

USA  Borl In 
Arm.  och  as) 

USASCTAF 

Arm.  och  (6) 

Allied  Command  Europe  (ACC) 

Aim.  OCH  (3) 

8th  USA.  Korpe  (14) 

ROK/US  Combined  Forces  Comnand  96301 
ATm.  CUSA-HHC-CFC/Cngr 

USA  Jspen  (USAAJ ) 

ATm*  AJEN-FE  96343 
ATm*  OCH- Honshu  96343 
ATm.  DEH-Ofc  fnava  96331 

Rocky  St.  Ares  80903 

Are*  Engineer.  AC  DO  Area  Office 
Arnold  Air  Force  Station,  m  37389 

Western  Area  Office.  CE 
Vandenberg  AFfl,  CA  93437 

416th  Engineer  Command  60623 
ATm.  Facilities  Engineer 

US  Military  Ac* dewy  10966 
ATTNt  Facilities  Engineer 
ATm  i  Dept  of  Geography  A 
Cowpvter  Science 
ATm*  OSCPCR/MACM-A 


mmc.  Arm.  o*ym-m  02172 

USA  ARROOK  61299 
ATm.  ORCIS-RI-I 
ATm.  OR  SAW- 1$ 

OAROON  -  01r..  Inst..  A  Svcs. 

ATTRi  OCH  123) 

OLA  ATm  1  OLA-WI  22314 

FORSOOK 

FORSOOK  Engineer.  ATm.  AFEM-FE 
ATm.  OCH  (23) 


MSI 

A I  IN  i  MSI  O-t  /H/n 
Aim.  labilities  1  n<|  lneer 
Mtrslnons  AMC  fUi'MO 
Walter  Reed  AMC  70012 

INSCON  -  Ch.  Inst).  01v. 

ATm.  Facilities  Engineer  (3) 

MOW 

Arm.  och  (3) 


Arm*  MTMC-SA  203  IS 

ATm.  Facilities  Engineer  (3) 

NAR ACCOM.  ATm.  DRONA-F  071160 

TARCOM,  Fac.  Dlv.  48090 

TRtfJOC 

HQ.  TRAOOC.  ATm.  ATEN-FE 
ATm.  DCH  (19) 

TSAROOM.  ATm.  STSAS-F  63120 


USACC 

ATm.  Facilities  Engineer  (2) 

WESTCOM 
AnN.  DCH 

Fort  Shatter  96856 
ATm.  AREN- IN 

SHAPE  09055 

ATm.  Survivability  Section.  CCS- OPS 
Infrastructure  Branch.  LANOA 

HQ  USEUCOM  09128 
Arm.  ecj  4/7-loc 

U.S.  Army,  Fort  Bolvoir  22060 
ATTN:  Canadian  Liaison  Officer 
ATTN:  Mater  Resources  Support  Center 
ATTN:  Engr  Studies  Center 
ATTN:  Engr  Topographic  Lab 
ATTN:  AT2A-0TE-SU 
ATTN:  AT2A-0TE-EM 

CRREL.  Arm.  Library  03755 

ETL,  ATm.  library  22060 

¥S,  ATm.  Library  39180 

HQ.  XVIII  Airborne  Corps  and 
Ft.  Bragg  28307 
ATm*  AFZA-FE-EE 

Chanvte  AFB.  IL  61068 
3345  CCS/OC.  Stop  27 

Norton  AFB  CA  92409 
ATm*  AFRCE-KX/OEE 

Tyndall  AFB.  FL  32403 

AFE SC/Engl neoring  4  Service  Lab 

NAFEC 

ATm*  R0T4E  Liaison  Offlco  (6) 

ATm.  Sr.  Tech.  FAC-03T  22332 
ATm.  Asst.  COR  RAO.  FAC-03  22332 

NCEl  93041 

Arm.  Library  (Code  L08A) 

Oefenae  Technical  Info.  Center  22314 
ATm.  OOA  (12) 

Engineering  Societies  Library 
Nev  York.  NY  10017 

National  Guard  Bureau  20310 
Installation  Division 

US  Government  Printing  Office  22304 
Receiving  Section/Depository  Copies  (2) 

US  Amy  Env.  Hygiene  Agency 

Arm.  msw-e  21010 

National  Bureau  of  Stoodor<i»  20760 


TOO 

1/9/84 


ENA  I tn  Distribution 


Chief  of  Engineers 
ATTN:  OALN-M'.C-i 
ATTN:  DAF.N-LCI  -  I* 

ATTN:  DAEN-ECE- I  12) 

ATTN:  DAEN-ZCF-il 
ATTN:  DAEN-CCZ-A 
ATTN:  DAEN-ZCE-D  (2) 

US  Army  Engineer  District 
New  York  10007 

ATTN:  Chief.  Design  Br 
Philadelphia  19106 
ATTN:  Chief.  NAPEN-E 
Baltimore  21203 
ATTN:  Chief.  Engr  Dlv 
Norfolk  Z3510 
ATTN:  Chief,  NAOEN-D 
Huntington  2572.1 
ATTN:  Chief.  fIRHED 
Wilmington  28401 
ATTN:  Chief,  SAWEN-0 
Savannah  31402 
ATTN:  Chief.  SASAS-L 
Mobile  36628 

ATTN:  Chief,  SAMEN-D 
Louisville  40201 

ATTN:  Chief,  Engr  Div 
St.  Paul  55 101 
ATTN:  Chief.  ED-D 
Chicago  60604 

ATTN.  Chief,  NCCPE-PES 
Rock  Island  61201 

ATTN:  Chief.  Engr  Oiv 
St.  Louis  63101 
ATTN:  Chief.  ED-D 
Omaha  6R102 
ATTN:  Chief,  Engr  Oiv 


New  Orleans  70160 
ATTN:  Chief,  LMNED-OG 
Little  Rock  72203 
ATTN:  Chief,  Engr  Div 
Tulsa  74102 
ATTN:  Chief,  Engr  Div 
Ft.  Worth  76102  (3) 

ATTN:  Chief.  SWFED-D 
San  Francisco  94105 
ATTN:  Chief,  Engr  Oiv 
Sacramento  95814 
ATTN:  Chief,  SPKED-D 
Far  East  96301 

ATTN:  Chief,  Engr  Oiv 
Seattle  98124 
ATTN:  Chief,  EN-D8-ST 
Walla  Walla  99362 

ATTN:  Chief,  Engr  Oiv 
Alaska  99501 
ATTN:  Chief,  NPASA-R 

US  Army  Engineer  Division 
New  England  02154 
ATTN:  Chief,  NEDED-T 
North  Atlantic  10007 
ATTN:  Chief,  NADEN-T 
Middle  East  (Rear)  22601 
ATTN:  Chief,  MEDEO-T 
South  Atlantic  30303 
ATTN:  Chief,  SAOEN-TS 
Huntsville  35807 
ATTN:  Chief,  HNOEO-CS 
ATTN:  Chief,  HNDED-SR 
Ohio  River  45201 
ATTN:  Chief,  Engr  Oiv 
Missouri  River  68101 
ATTN:  Chief,  MROED-T 
Southwestern  75202 
ATTN:  Chief,  SWOEO-T 
South  Pacific  94111 
ATTN:  Chief,  SPOEO-TG 
Pacific  Ocean  96858 
ATTN:  Chief,  Engr  Oiv 
North  Pacific  97208 


ij  tn  US  Army  94129  N 

AT1N  AFKf.-IH 

N 

7th  Army  Combined  Arms  Trrig.  Cntr.  U94U7 
ATTN-  AETTM-HRO-FHO  V 


•  .  mm\n  • 


Armament  1  Oev.  Command  21005 
ATTN:  0R9AR-BLT 

USA  ARRAOCOM  07801 
4TrN ■  ORDAR-LCA-OK 

DARCOM  22333 
ATTN-  DRCPA-E 
ATTN  ORCIS-A 

TRAOOC 

Ft.  Monroe.  VA  23651 

Ft.  Clayton,  Canal  Zone  34004 
ATTN:  OF  AE 

Ft.  Detrick,  MO  21701 

Ft.  Leavenwortn,  KS  66027 
ATTN:  ATZLCA-5A 

Ft.  McPherson,  GA  30330  (2) 

Ft.  Monroe,  VA  23651  (6) 

Ft.  Rucker,  AL  36360  (2) 

Aberdeen  Proving  Ground,  MO  21005 
ATTN:  DROAR-BLL 
ATTN:  STEAP-MT-E 

Human  Engineering  Lab.  21005  (2) 

USA-WES  39181 

Army  Environmental  Hygiene  Agency  21005 

Naval  Air  Station  92135 
ATTN:  Code  661 


NAVFAC  22332  (2) 

Naval  Air  Systems  Coraand  20360 

US  Naval  Oceanographic  Office  39522 

Naval  Surface  Weapons  Center  22485 
ATTN:  N-43 

Naval  Undersea  Center,  Code  401  92152  (2) 

Bolling  AFB,  DC  20332 
AF/LEEEU 

Patrick  AFB,  FL  32925 
ATTN :  XRQ 

Tyndall  AFB,  FL  32403 
AFESC/TST 

Wright-Patterson  AFB,  OH  45433  (3) 

Building  Research  Advisory  Board  20418 

Transportation  Research  Board  20418 

Oept  of  Housing  and  Urban  Development  20410 

Dept  of  Transportation  Library  20590 

Illinois  EPA  62706  (2) 

Federal  Aviation  AAiinlstration  20591 

Federal  Highway  Administration  22201 
Region  15 


NASA  23365  (2) 

National  Bureau  of  Standards  2U234 

7 

Office  of  Noise  Abatement  20590 
ATTN:  Office  of  Secretary 

USA  logistics  Management  Center  23801 

Airports  and  Construction  Services  Dir 
Ottawa,  Ontario,  Canada  IC1A  0N8 

Division  of  Building  Research 
Ottawa,  Ontario,  Canada  K1A  0R6 

National  Defense  HQDA 
Ottawa,  Ontario,  Canada  K1A  0K2 


>.*  -.1 

Xv- v 1 


Schomer,  Paul  D. 

Strategies  for  and  validity  of  noise  monitoring  in  the  vicinity  of 
civilian  airfields  and  Army  installations  /  by  Paul  D.  Schomer,  Richard  E 
DeVor.  —  Champaign,  Ill  :  Construction  Engineering  Research  Laboratory, 
1983. 

48  p,  (Technical  report  /  Construction  Engineering  Research  Laboratory 
N-166) 

1.  Airport  noise — measurement.  1.  DeVor,  Richard  E.  II.  Title. 
III.  Series;  Technical  report  (Construction  Engineering  Research  labor  a 
tory )  ;  N-166. 


TO 


